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A study of the D*j(2317) + and D s j (2460) + mesons in inclusive cc production is presented using 
232 fb _1 of data collected by the BABAR experiment near yfs = 10.6 GeV. Final states consisting 
of a Df meson along with one or more 7r°, it , or 7 particles are considered. Estimates of the mass 



and limits on the width are provided for both mesons and for the D s i(2536) + meson, 
also performed for neutral and doubly-charged partners of the D*j (2317) + meson. 

PACS numbers: 14.40.Lb, 13.25.Ft, 13.20.Fc 
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I. INTRODUCTION 



The D*j(2317) + meson, discovered by this collabora- 
tion [1] and confirmed by others [2, 3], does not conform 
to conventional models of cs meson spectroscopy. In- 
cluded with this discovery were suggestions of a second 
state, the D s j (2460) + meson. This meson, observed by 
CLEO [2] and confirmed by this collaboration [4] and 
BELLE [3], has a mass that is also lower than expecta- 
tions [5-8] . Because the masses of these two states are 
so unusual there has been speculation [9, 10] that both 
possess an exotic, four-quark component. The possibil- 
ity that the D* sJ {2317)+ and D sJ (2460)+ are exotic has 
attracted considerable experimental and theoretical in- 
terest and has focused renewed attention on the subject 
of charmed-meson spectroscopy in general. 

Presented in this paper is an updated analysis of these 
two states using 232 fb~ of e + e~ — > cc data collected 
by the BABAR experiment. From this analysis new esti- 
mates of the L>* J (2317)+ and D sJ (2460) + masses, limits 
on their intrinsic widths, calculations on their production 
cross sections, and the branching ratios of .D s j(2460) + 
decays to -D+7 and D+7r + 7r~ with respect to its decay 
to D+7r°7 are presented. 

These measurements are performed by fitting the in- 
variant mass spectra of combinations [11] of D^~tt°, -D+7, 
£)+7r°7, -D+7r 7r°, Z?+77, and Dfir + 7r~ particles. Com- 
binations of D+tt + and D+tt~ are also studied to search 
for new states. The mass spectrum of each final-state 
combination is studied in detail. In particular, features 
in the spectra that arise from reflections of other cs meson 
decays are individually identified and modeled. The anal- 
ysis of the Z?+7r°7 final state includes an explicit search 
for the two most likely sub-resonant decay channels for 
.D s j(2460) + meson decay. 

This paper is organized as follows. First, the current 
status of the D* sJ (2317)+ and L> sJ (2460) + 
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mesons is re- 
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didates is then described, including an estimate of Dj 
yield in terms of the Df — > cf>TT + branching fraction. 
Each combination of final-state particle species is then 
discussed individually. The paper finishes with a sum- 
mary of results and conclusions. 
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FIG. 1: The cs meson spectrum, as predicted by Godfrey and 
Isgur [5] (solid lines) and Di Pierro and Eichten [8] (dashed 
lines) and as observed by experiment (points). The DK and 
D* K mass thresholds are indicated by the horizontal lines 
spanning the width of the plot. 



II. REVIEW OF THE D* J (2317)+ AND 

£> sJ (2460) + 



Much of the theoretical work on the cs system has been 
performed in the limit of heavy c quark mass using po- 
tential models [5-8] that treat the cs system much like a 
hydrogen atom. Prior to the discovery of the D*j(2317) + 
meson, such models were successful at explaining the 
masses of all known D and D s states and even predicting, 
to good accuracy, the masses of many D mesons (includ- 
ing the D s i(2536) + and D S 2(2573) + ) before they were 
observed (see Fig. 1). Several of the predicted D s states 
were not confirmed experimentally, notably the lowest 
mass J p = + state (at around 2.48 GeV/c 2 ) and the sec- 
ond lowest mass J p = 1 + state (at around 2.58 GeV/c 2 ). 
Since the predicted widths of these two states were large, 
they would be hard to observe, and thus the lack of ex- 
perimental evidence was not a concern. 

The D* 7 (2317) + meson has been observed in the decay 
L»* 7 (2317)+ -» D+tt° [1-3, 12, 13]. The mass is mea- 
sured to be around 2.32 GeV/c 2 , which is below the DK 
threshold. Thus, this particle is forced to decay either 
electromagnetically, of which there is no experimental ev- 
idence, or through the observed isospin-violating Dfir 
strong decay. The intrinsic width is small enough that 
only upper limits have been measured (the best limit 
previous to this paper being T < 4.6 MeV at 95% CL 
as established by BELLE [3]). If the D* 7 (2317) + is the 



missing + cs meson state, the narrow width could be ex- 
plained by the lack of an isospin-conserving strong decay 
channel. The low mass (160 MeV/c 2 below expectations) 
is more surprising and has led to the speculation that the 
D*j(2317) + does not belong to the D+ meson family at 
all but is instead some type of exotic particle, such as a 
four-quark state [9] . 

The D s j(2460) + meson has been observed decaying 
to D+7r°7 [2-4, 12, 13], D+TT+n- [3], and B>+ 7 [3, 12, 
13]. The intrinsic width is small enough that only upper 
limits have been measured (the best limit previous to this 
paper being T < 5.5 MeV at 95% CL as established by 
BELLE [3]). The D+7 decay implies a spin of at least 
one, and so it is natural to assume that the D s j(2'±60) + is 
the missing 1 + cs meson state. Like the D*j(2317) + , the 
D s j(24Q0) + is substantially lower in mass than predicted 
for the normal cs meson. This suggests that a similar 
mechanism is deflating the masses of both mesons, or 
that both the states belong to the same family of exotic 
particles. 

The spin-parity of the B>* 7 (2317)+ and D sJ (2460) + 
mesons has not been firmly established. The decay mode 
of the Z?*j (2317) + alone implies a spin-parity assign- 
ment from the natural J p series {0 + , 1~, 2 + , . . . }, as- 
suming parity conservation. Because of the low mass, 
the assignment J p = + seems most reasonable, al- 
though experimental data have not ruled out higher 
spin. It is not clear whether electromagnetic decays such 
as .D*(2112) + 7 can compete with the strong decay to 
-D+7T , even with isospin violation. Thus, the absence 
of experimental evidence for radiative decays such as 
D* 7 (2317) + -> D*(2112)+7 is not conclusive. 

Experimental evidence for the spin-parity of the 
.D s j(2460) + meson is somewhat stronger. The observa- 
tion of the decay to Dfj alone rules out J = 0. Decay 
distribution studies in B ^ D sJ (2460) + D { s * ) ^ [12, 13] 
favor the assignment J = 1. Decays to either Dfn°, 
D°K + , or D + K° would be favored if they were al- 
lowed. Since these decay channels are not observed, this 
suggests, when combined with the other observations, 
the assignment J p = 1 + . In this case, the decay to 
D*j(2317) + 7 is allowed, but it may be small in compar- 
ison to the -D+7 decay mode. 

Table I lists various possible decay channels for the 
D* 7 (2317)+ and D Si7 (2460) + mesons. Several of these 
decays are forbidden assuming the spin-parity assign- 
ments discussed above. 



III. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were recorded with the 
BABAR detector at the PEP-II asymmetric-energy stor- 
age rings and correspond to an integrated luminosity of 
232 ftV 1 collected on or just below the T(45) resonance. 

A detailed description of the BABAR detector is 
presented elsewhere [14]. Charged particles are de- 
tected with a five-layer, double-sided silicon vertex 



TABLE I: A list of various decay channels and whether they 
have been seen, are allowed, or are forbidden in the decay 
of the D* sJ (2317) + and D sJ (2460) + mesons. The predictions 
assume a spin-parity assignment of J p — + and 1 + , respec- 
tively. 



Decay Channel 


£> s *j(2317)+ 


D sJ (2460)+ 


Dt^ 


Seen 


Forbidden 


D+7 


Forbidden 


Seen 


D+7r°7 (a) 


Allowed 


Allowed 


D:(2112) + 7t° 


Forbidden 


Seen 


D:. 7 (2317) + 7 


— 


Allowed 


D+ttV 


Forbidden 


Allowed 


D+77 (a) 


Allowed 


Allowed 


D s *(2112)+7 


Allowed 


Allowed 


D+TT + TT" 


Forbidden 


Seen 



(a) Non-resonant only 



tracker (SVT) and a 40-layer drift chamber (DCH) us- 
ing a helium-isobutane gas mixture, placed in a 1.5-T 
solenoidal field produced by a superconducting magnet. 
The charged-particle momentum resolution is approxi- 
mately (8 P t/pt) 2 = (0.0013 p T ) 2 + (0.0045) 2 , where p T 
is the transverse momentum in GcV/c. The SVT, with 
a typical single- hit resolution of 10 /im, measures the im- 
pact parameters of charged-particle tracks in both the 
plane transverse to the beam direction and along the 
beam. Charged-particle types are identified from the 
ionization energy loss (dE/dx) measured in the DCH 
and SVT, and from the Cherenkov radiation detected 
in a ring- imaging Cherenkov device (DIRC). Photons 
are detected by a CsI(Tl) electromagnetic calorimeter 
(EMC) with an energy resolution <j(E)/E = 0.023 • 
(E/ GcV)~ 1/4 ® 0.019. The return yoke of the supercon- 
ducting coil is instrumented with resistive plate chambers 
(IFR) for the identification of muons and the detection 
of neutral hadrons. 



IV. CANDIDATE RECONSTRUCTION 



The goal of this analysis is to study the possible decay 
of the D*j(2317)+ and D sJ (2460)+ mesons into the fi- 
nal states listed in Table I. These decay channels consist 
of one Df meson combined with up to two additional 
particles selected from tt°, it , and 7. The first step 
in this analysis is to identify D+ mesons in the BABAR 
data. For each resulting Df candidate a search is per- 
formed for associated 7r°, 7, and ir particles. Signals 
from D* 7 (2317)+ and £> sJ (2460) + decay are isolated us- 
ing the invariant mass of the desired combination of par- 
ticle species. 

The Df — ► K + K~n + decay mode is used to select a 
high-statistics sample of Df meson candidates. Each K^ 
and -k^- candidate is separated from other charged par- 
ticle species by a likelihood-based particle identification 
algorithm based on the Cherenkov-photon information 
from the DIRC together with dE/dx measurements from 
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FIG. 2: (a) The invariant mass spectrum of Df candidates 
before (top histogram) and after (bottom points) applying 
subresonant (f>ix + and K*K + selection. The light (dark) ar- 
eas indicate the signal (sideband) regions, (b) The invariant 
77 mass of ir° candidates before (top histogram) and after 
(bottom points) applying the tv° veto described in the text. 
The light histogram indicates those candidates that pass the 
X 2 requirement. The curve in (a) is the \ 2 fit described in 
the text. 



the SVT and DCH. A geometrical fit to a common vertex 
is applied to each K + K~ir + combination. An acceptable 
K + K~tt + candidate must have a fit probability greater 
than 0.1% and a trajectory consistent with originating 
from the e + e~ luminous region. To reduce combinato- 
rial background, each K + K~tt + candidate must have a 
momentum p* in the e + e~ ccntcr-of-mass frame greater 
than 2.2 GeV/c 2 , a requirement that also removes nearly 
all contributions from B-meson decay. Background from 
D° — > K + K~, which is evident from the corresponding 
K + K~ mass distribution, is removed by requiring that 
the K + K~ mass be less than 1.84 GeV/c 2 . 

The upper histogram in Fig. 2(a) shows the K + K~ir + 
mass distribution for all candidates. A clear Df signal 
is seen. To reduce the background further, only those 
candidates with K + K~ mass within 10 McV/c 2 of the 
0(1020) mass or with K~n + mass within 50 MeV/c 2 
of the if* (892) mass arc retained; these densely pop- 
ulated regions in the Df Dalitz plot do not overlap 
(see Fig. 3). The decay products of the vector particles 
0(1020) and J?* (892) exhibit the expected cos 2 9 h be- 
havior required by conservation of angular momentum, 
where Oh is the helicity angle. The signal-to-background 
ratio is further improved by requiring |cos^^| > 0.5. 
The lower histogram of Fig. 2(a) shows the net effect 
of these additional selection criteria. The Z)+ signal 
(1.954 < m(K+K-n+) < 1.981 GeV/c 2 ) and sideband 
(1.912 < m(K+K-TT+) < 1.934 GeV/c 2 and 1.998 < 
m(K+K-n + ) < 2.020 GeV/c 2 ) regions are shaded. This 
distribution can be reasonably modeled in a \ 2 fit by the 
sum of two Gaussian distributions with a common mean 
(hereafter referred to as a double Gaussian) on top of a 
quadratic background. The result of this fit is a Df sig- 
nal peak consisting of approximately 410 000 decays and 
a mass of 1967.8 MeV/c 2 with negligible statistical error. 
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FIG. 3: An illustration of the Dt sub-resonant selection re- 
quirements. The light shaded area is the kinematic range 
for Df — > K + K~tt + decay. The vertical line with an 
arrow represents the selection requirement used to remove 
D° — * K + K~ decay. The four dark regions indicate those 
portions of the Df phase space used for final candidate selec- 
tion, corresponding to Df — > 4>ir + and Df — » K* K + decay. 
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FIG. 4: The sideband subtracted p* spectrum for Df candi- 
dates after all selection requirements are met. 



The approximate p* distribution for selected D+ 
mesons can be obtained by simple sideband subtraction, 
assuming linear background behavior under the Df . The 
result is shown in Fig. 4. 

The final list of D+ candidates are those that lie within 
the signal window. For each such candidate, the mo- 
mentum vector is calculated from the simple addition 
of K + , K~, and ir + momentum vectors. The energy 
is chosen to reproduce the PDG value for the D+ mass 
(1968.5 ±0.5) MeV/c 2 [15]. 

It is assumed that the D*j(2317) + and £> sJ (2460) + 
mesons have lifetimes that are too small to be resolved by 
the detector. Thus, the most likely point of Z?*j(2317) + 
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FIG. 5: The D^"-sideband subtracted p* spectrum for associ- 
ated (a) ty , (b) 7r , and (c) 7, after the selection requirements 
described in the text are fulfilled. 



and .D s j(2460) + decay for each Df candidate is chosen 
to be the interaction point (IP), calculated from the in- 
tersection of the trajectory of the candidate and the e + e~ 
luminous region. To produce a list of ir particles that 
could arise from D* sJ (2317) + or £> sJ (2460) + decay, the 
trajectories of all ^ candidates that are not daughters 
of the D+ candidate are constrained to the IP using a 
geometric vertex fit. The approximate p* spectrum of 
these candidates associated with real Df mesons can be 
obtained by using simple Df sideband subtraction. The 
result is shown in Fig. 5(a). 

The selection of 7 and 7r° candidates is a two-step pro- 
cess. The first step is the selection of a fiducial list of 
7 and 7T° candidates. The fiducial list of 7 candidates is 
constructed from energy clusters in the EMC with ener- 
gies above 100 MeV and not associated with a charged 
track. The energy centroid in the EMC combined with 
the IP position is used to calculate the 7 momentum di- 
rection. Each fiducial ir° candidate is constructed from 
a pair of 7 particles in the 7 fiducial list. This 7 pair is 
combined using a kinematic fit assuming a tt° mass. The 
resulting ir° momentum is required to be greater than 
150 McV/c. The 77 invariant mass spectrum from this 
7t° selection is shown in the top histogram of Fig. 2(b). 
To produce the final fiducial list of ir° candidates, the x 2 
probability of the kinematic fit is required to be greater 
than 2%. 

The final list of 7 candidates consists of any 7 in the 
fiducial list that is not used in the construction of any 
7r° in the w° fiducial list. The final list of ir° candidates 
consists of any w° in the fiducial list that does not share 
a 7 with any other ir° in the fiducial list. The 77 invari- 



ant mass distribution of the final list of tt° candidates is 
shown in the bottom histograms of Fig. 2(b), both be- 
fore and after applying the \ 2 probability requirement. 
The approximate p* spectrum of the final list of 7 and it 
candidates as determined using Df sideband subtraction 
is shown in Figs. 5(b) and 5(c). 



V. MONTE CARLO SIMULATION 

Monte Carlo (MC) simulation is used for the following 
purposes in this paper: 

• To calculate signal efficiencies. 

• To provide independent estimates of background 
levels. 

• To characterize the reconstructed mass distribution 
of the signal. 

• To predict the behavior of various specific types of 
backgrounds (commonly referred to as reflections) 
produced when the mass distribution from an es- 
tablished decay mode is distorted by the loss of one 
or more final-state particles or the addition of one 
or more unassociated particles. 

Various sets of MC events were generated. For the pur- 
poses of understanding signal efficiencies, signal shapes, 
and reflections, individual MC sets of 500 000 decays were 
generated for each known decay mode of the D*j(23l7) + 
and D sJ (2460)+ mesons. In addition, MC sets of 250 000 
decays were produced for each hypothetical D* sJ (2317)+ 
and _D s j(2460) + decay as needed. Finally, a set of 
e + e~ — » cc events, corresponding to an integrated lumi- 
nosity of approximately 80 fb~ , was generated to study 
sources of combinatorial background. 

Each MC set was processed by the same reconstruction 
and selection algorithms used for the data. Independent 
tests of the detector simulation have demonstrated an 
accurate reproduction of charged particle detection effi- 
ciency. The systematic uncertainty from these tests is 
estimated to be 1.3% for each charged track. Since the 
decay D+ — > K + K~ir + involves three charged tracks, 
the systematic uncertainty in Df efficiency from the sim- 
ulation alone is estimated to be 3.9%. 

The simulation of D+ — > K + K~w + decay was de- 
signed to match approximately the known Dalitz struc- 
ture. MC events were reweightcd to match more precisely 
the relative <fiir + and K*K + yields observed in the data. 

The simulation assumes a D*j(23f7) + and 
D s ,7(2460)+ intrinsic width of V = 0.1 MeV. All 



_D*j(2317) + and D s j(2460) + final states are generated 
using phase space. The generated p* distribution of 
-0^(2317)+ and D sJ (2460)+ mesons in the MC simula- 
tion was adjusted to roughly reproduce observations. 



VI. ABSOLUTE D+ YIELD 



In order to calculate D*j(23l7) + and D sJ (2A60)+ pro- 
duction cross sections, it is necessary to provide an esti- 
mate of absolute D+ selection efficiency. Since this anal- 
ysis uses Df — * K + K~n + decay, the approach is to nor- 
malize Df yield with respect to the Df — > (f>ir + , <j> — > 
K + K~ branching fraction, the world average of which is 
1.8 ± 0.4% [15]. To perform this normalization correctly, 
the following must be accounted for: 

• The K*K + portion of the Df sample. 

• Non-resonant K + K ~tt + background under the <j) 
peak. 

• The fraction of the (f> signal that falls outside of the 
K + K~ mass selection and Oh requirements. 

The K*K + selection represents approximately 48% of 
the total D+ sample. An inspection of the p* distribu- 
tion of the <fiir + and K*K + subsamples indicates that 
this fraction is, to a good approximation, independent of 
p* . Therefore, a constant factor is sufficient to account 
for the contribution from the K*K + portions of the D+ 
sample. 

Shown in Fig. 6 is the Dj~-sideband subtracted K + K~ 
invariant mass spectrum for all Df candidates before ap- 



and K*K + selection requirements. 

,2 



A 



plying the </> 

prominent (f> peak is observed. A binned \ 2 fit to this 
spectrum is used to extract both the fraction of (f>ir + de- 
cays that fall outside the <f> selection window and the 
number of non-0 decays that leak inside. This fit is de- 
scribed below. 

To model the signal portion of the K + K~ mass spec- 
trum, the 4> - * K + K~ line shape a(m) can be reasonably 
well described (ignoring potential interference effects) by 
a relativistic Breit-Wigncr function: 



a(m) = 



toF(to) 



(m 2 -ml) +mgr t 2 ot (m) 



(1) 



where m = (1019.456 ± 0.020) MeV/c 2 is the intrinsic 
4> mass [15]. The mass-dependent width T{m) can be 
approximated by: 



T(m) = 0.493r 



m 



x 3 1 



9o 



2d2 



q$R 



1 + q 2 R 2 ' 



(2) 



where To = 4.26 MeV is the intrinsic width, 0.493 is the 
branching fraction for this decay mode, R = 3 GeV~ is 
an effective <j) radius (to control the tails), and q (qo) is 
the total three-momentum of the K^ decay products in 
the (f> center-of-mass frame assuming an effective (f> mass 
of m (mo): 



q = \lm 2 /A — m 2 K and qo = \ TOq/4 — m 2 K . 



(3) 



A reasonable approximation of the total width in- 
cludes the three dominant decay modes (K + K~, KsKl, 
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FIG. 6: The DJ -sideband subtracted K K~ invariant mass 
spectrum near the (f> mass for the Df sample obtained before 
applying the <fm + and K* K + selection requirements. The 
curve is the fit described in the text. The dashed line is the 
portion of the fit attributed to contributions from other than (j> 
decay. The vertical lines indicate the 4> mass selection window. 



7T + vT 7T°): 



r t ot(m) = T(m) + r'(m) + 0.155F 



(4) 



where r'(m) is calculated in the same manner as T(m) 
but using the Ks/Kl mass and branching fraction 
(0.337) and the width for the three-pion decay mode 
(which is well above threshold) is treated as a constant. 
The fit function P{m) used to describe the K + K~ 
mass spectrum of Fig. 6 is: 



P(m) = S^m) + C{m) 



(5) 



where S^,{m) is the (f> signal shape and C(m) represents 
non-resonant contributions. The empirical form used for 
C(m) is a four-parameter threshold function: 



C{m) = 



to < a\ 

C'{m) to > a\ 

C"(to) = (l + TO- ai ) a2 (l-cxp(-(TO-a )/a3)) 



(6) 



The form used for S^m) is the Breit-Wigncr function of 
Eq. 1 smeared by a Gaussian: 



S^(to) 



N 



SV2ir 



<j(m')e-( m - m '¥ /2s2 dm' , (7) 



where N is an overall normalization. In the fit to this 
function, the value of To is kept fixed to the PDG value 
but Too and 5 are allowed to vary. 



Despite its simplicity, the function of Eq. 5 describes 
the K + K~ spectrum quite well, as shown in Fig. 6. 
The value of too produced by the fit is slightly lower 
[(—56 ± 6) keV/c 2 , statistical error only] than the PDG 
average [15]. To determine a correction factor for the D + 
yield, the total <fi yield (calculated from the integral of the 
signal line shape determined by the fit up to a K + K~ 
mass of 1.1 GeV/c 2 ) can be compared to the total num- 
ber of candidates which fall inside the cf> mass window. 
The result is a correction factor of 1.09, with negligible 
statistical uncertainty. 

To test the above calculation, the fit is repeated on a 
D + sample that includes the |cos6^| > 0.5 requirement 
discussed earlier for the <\>. The change in cf> integrated 
yield is consistent with a cos 2 9h distribution. The mea- 
sured correction factor increases to 1.10. The difference 
between this value and 1.09 is taken as a systematic un- 
certainty. 

Other systematic checks performed include increasing 
the range in K + K~ mass of the (f> line shape integration 
and changing the value of R from 1 to 5 GeV/c 2 . The to- 
tal uncertainty in the 1.09 correction factor is found to be 
a 0.043 (3.9% relative), as calculated from a quadrature 
sum. 



VII. CANDIDATE SELECTION OPTIMIZATION 

This paper explores the eight final-state combinations 
shown in Table II, each involving a D + meson and up to 
two total of TT+ , 7T°, and/or 7 particles. For each combi- 
nation it is necessary to distinguish possible signals from 
D* 7 (2317)+ and L> sJ (2460)+ decay from combinatorial 
background. The separation of signal and background 
is made more distinct if additional candidate selection 
requirements are imposed. This section discusses those 
additional requirements. 

In four of the final-state combinations (D + tt°, D + tt°j, 
Dfj, and D+7r + 7r~) a signal is expected. An estimate 
of signal significance is calculated in these cases based on 
expected signal and background rates, the former calcu- 
lated from previously published branching ratio measure- 
ments combined with the appropriate MC sample. For 
the remainder of the final states, an estimate of signal 
sensitivity is calculated for the hypothetical D* sJ (2317) + 
and D s j(2460) + meson decay. This sensitivity calcula- 
tion is based on signal efficiency, determined using MC 
samples, and expected background levels. 

To avoid potential biases, both the signal significance 
and sensitivity estimates are calculated solely using MC 
samples. 

The following selection requirements are adjusted in 
order to produce optimal values of signal significance and 
sensitivity: 

• A minimum total center-of-mass momentum p* . 

• A minimum energy (for 7) and/or momentum cal- 
culated in the laboratory (for ir°) or center-of-mass 



TABLE II: Selection requirements for the final states studied 
in this paper, the resulting number of events, and the ap- 
proximate efficiency for a D s j (2460) + signal. The selection 
requirements are specified either in the laboratory (Lab) or 
center-of-mass (CMS) coordinate systems. 





Minimum Requirements 








7 Energy 


7r u Mom. 


7T* Mom. 






Lab 


Lab 


CMS 


Sample 


Effic. 


Final State 


(MeV) 


( MeV/c) 


( MeV/c) 


Size 


(%) 


Dt^ 


— 


350 


— 


87 320 


6.4 


D+~i 


500 


— 


— 


133 398 


12 


D+7T°7 


135 


400 


— 


170 341 


2.4 


D + 7V°7T° 


— 


250 


— 


17 437 


0.4 


D+jj 


170 


— 


— 


575 765 


7.9 


D+K+ 


— 


— 


300 


143 149 


13 


D + 7Y- 


— 


— 


300 


219 466 


13 


D+TY+n- 


— 


— 


250 


154 496 


6.8 



(for 7r ) frame of reference. 

For the minimum p* , it is important to choose the same 
value for all final-state combinations in order to mini- 
mize systematic uncertainties in the branching ratios. A 
minimum value of p* > 3.2 GeV/c is chosen as a rea- 
sonable compromise. Values for the remaining selection 
requirements are chosen separately for each final-state 
combination. The results are listed in Table II. 

The MC samples can be used to estimate the approx- 
imate efficiency for detecting a signal with a p* of at 
least 3.2 GeV/c after applying the above selection re- 
quirements. The resulting efficiencies vary between 0.4% 
and 13% (see Table II). 



VIII. CROSS SECTION NOTATION 

To report production yields of a particular cs meson 
Dy to a particular final state D + X, the following quan- 
tity a(D Y -> D+X) is defined: 



a(D Y -► D+X) = a{e+e 
x B{D Y 
x B(D+ 



-* D Y ,p Y >3.2 GeV/c) 

DfX) 
cjm+^^K+K-), (8) 



where the cross section a is defined for a center-of-mass 
momentum p* above 3.2 GeV/c. The quantity a is cal- 
culated by taking the number of Dy —* D + X decays 
observed in the data, correcting for efficiency using the 
appropriate MC sample (restricted to p* > 3.2 GeV/c), 
correcting for the relative D + — > (f)n + yield as calculated 
in Section VI, and dividing by the luminosity (232 fb _ ). 
A relative systematic uncertainty of 1.2% is introduced 
to account for the uncertainty in the absolute luminosity. 
There is no attempt to correct the cross section for ra- 
diative effects (such as initial-state radiation). Since a 
reasonably accurate representation of such radiative ef- 
fects is included in our MC samples, the calculation of se- 



lection efficiencies from these samples is accurate enough 
for the purposes of this paper. 
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IX. THE D+tt° FINAL STATE 

Shown in Fig. 7 is the invariant mass distribution of 
the -Dj~7r° combinations after all selection requirements 
are fulfilled. Signals from L>*(2112)+ and D* sJ (2317) + 
decay arc evident. An unbinncd likelihood fit is applied 
to this mass distribution in order to extract the parame- 
ters and yield of the D*j (2317) + signal and upper limits 
on Z? s ,/(2460) + decay. The likelihood fit includes six dis- 
tinct sources of Dfir combinations: 

• D* sJ (2317)+ -► D+ir" decay 

• D sJ (2460)+ -> D+ir° decay (hypothetical). 

• D*(2112)+ -► D+7T decay. 

• A reflection from D*(2112) + — » Dfj decay in 
which an unassociated 7 particle is added to form 
a false ir° candidate. 

• A reflection from D sJ (2460)+ -► D*(2112)+tt° de- 
cay in which the 7 from the D*(2112) + decay is 
missing. 

• Combinatorial background from unassociated D^ 
and 7T° mesons. 

The probability density function (PDF) used to describe 
the mass distribution of each of these sources is described 
below. 

As shown in Fig. 8(a), the reconstructed mass distri- 
bution of the L>* 7 (2317)+ -► L>+tt° decay, as predicted 
by MC, has non-Gaussian tails and is slightly asymmet- 
ric. To describe this shape, the MC sample is fit to a 
modified Lorentzian function F^{m): 



F L {m) 



03- 



ll + a 4 <5 + a 5 (5 3 



(1 + S 2 ) ae 
S = (m — 01) /02 , 



(9) 



where a\ and a 2 correspond roughly to a mean and width, 
respectively. This function is simply a convenient pa- 
rameterization of detector resolution. The fit results are 
shown in Fig. 8(a). A similar procedure is used for the 
hypothetical L> sJ (2460)+ -> L>+tt° decay (Fig. 8(b)). 

The 350 MeV/c tt° momentum requirement removes 
the majority of L>*(2112)+ -> L>+tt decays. The re- 
maining signal is modeled using a distribution J(m) of 
Gaussians: 



J(m) = a-[ 



1 



a^a- 



-exp [-(to - a 2 f/2a 2 } da . (10) 



The parameters a 3 and 04 of this function are determined 
using a fit to a suitable MC sample. The mean mass 02 
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FIG. 7: The invariant mass distribution for (solid points) 
Df-K candidates and (open points) the equivalent using the 
Df sidebands. The curve represents the likelihood fit de- 
scribed in the text. Included in this fit is (light shade) a con- 
tribution from combinatorial background and (dark shade) 
the reflection from D sJ (2460)+ -> D*(2112) + 7r° decay. The 
insert highlights the details near the D*j(2317) + mass. 



is set equal to 2112.9 MeV/c 2 (0.5 MeV/c 2 higher than 
the PDG value [15]) to match the data. 

A reflection in D+tt produced by D*(2112)+ -► D+7 
decay appears as a broad distribution peaking at a mass 
of approximately 2.17 GeV/c 2 . This reflection is pro- 
duced by fake ir° candidates consisting of the 7 parti- 
cle from £>*(2112) + decay combined with unassociated 
7 candidates. Kinematics limit this reflection to Dfir 
masses above the quadrature sum of Df and ir° meson 
masses (approximately 2.1167 GeV/c 2 ). This distribu- 
tion falls gradually as mass is increased due to the rapidly 
falling inclusive 7 energy spectrum. Detector resolution 
tends to smear the lower kinematic mass limit. To model 
this reflection, a quadratic function with a sharp lower 
mass cut off is convoluted with a Gaussian distribution. 
The parameters of this function are determined directly 
from the D+ir° data sample. 

The D s j(2460) + reflection requires careful attention 
because it appears directly under the D* sJ (231 7) + sig- 
nal. This reflection is produced by the Dfir projection 
of L> sJ (2460)+ -► D*(2112)+7r° decay (in which the 7 
from D*(2112) + decay is ignored). A kinematic calcu- 
lation (Fig. 8(c)) of the £> sJ (2460) + Dalitz distribution 
predicts that this reflection, at the limit of perfect resolu- 
tion and efficiency, is a flat distribution in mass squared 
centered at a Dfn mass of 2313.4 MeV/c 2 with a full 
width of 41.3 MeV/c 2 (assuming a L» sJ (2460) + mass of 
2458.0 MeV/c 2 ). 



11 




800 -(b) 



600 



200 



2.3 2.35 2.4 
0; tt° Mass (GeV/c 2 ) 



2.-1 



f 



2.45 



" 100 



2.5 



Ds tt° Mass (GeV/c 2 



2.4 



2.2 



-(<0 




D;(21 12) 


"ill, 




I 



400 


""" ' » 




300 


" f 


\ 


200 


f 

- 1 


\ 


100 


\ 




; ti 

r I ,~jf , I 


• 
,\ 



2.2 2.4 

D>° Mass (GeV/c 2 ) 



2.25 



2.3 2.35 
D s * tt° Mass (GeV/c 2 ) 



FIG. 8: The reconstructed Dfir invariant mass spectrum 
from (a) D* sJ (2317) + and (b) D sJ (2460) + MC samples. The 
fit function of Eq. 9 is overlaid. (c) The projection of 
-D+7T mass for subresonant D s j(2A60) + decay through the 
D*(2112) + meson is restricted to a narrow range centered 
around 2313.4 MeV/c 2 . (d) The reconstructed Dfir invari- 
ant mass spectrum for the D s j(2460) + reflection from MC 
simulation. The solid curve is the fit function. The dashed 
curve is the same fit function with the Gaussian smearing 
removed. 



The D s ,/(2460) + reflection is flat in mass squared only 
if the 7T° efficiency is constant. In practice this is not 
the case, as illustrated by MC simulation (Fig. 8(d)). 
To accommodate the non-constant efficiency, the D+ir° 
mass distribution from the MC sample is fit to a function 
consisting of a bounded quadratic function smeared by a 
double Gaussian. The result of the fit is shown in Fig. 8d. 

The threshold function C(m) of Eq. 6 is used to repre- 
sent the mass spectrum from combinatorial background 
where the threshold value a\ is fixed to 2103.5 MeV/c 2 , 
the sum of the assumed Df and it masses. The remain- 
ing parameters of C(m) are determined directly from the 
data. 

The results of the likelihood fit to the D+ir° mass spec- 
trum is shown in Fig. 7. In this fit, the size, shape, and 
mean mass of the D s j(2460) + reflection are fixed to val- 
ues consistent with the yield and mass results determined 
in Section XI of this paper. The yield of Z?*(2112) + and 
D* sJ (2317) + decay and the D*j (2317) + mass is allowed 
to vary to best match the data. A D*j(2317) + mass of 
(2319. 6±0. 2) MeV/c 2 is obtained (statistical error only) . 
A total of 3180 ± 80 L»* / (2317)+ decays are found. 

The fit includes a hypothetical contribution from 
.D s j(2460) + — > Df 7T° in the form of a line shape of fixed 
shape and mass. The result is a yield of — 40 ± 50 (statis- 
tical errors only). The size of this yield is small enough 
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FIG. 9: Fit results near the D*j(2317)+ peak for the D+ir° 
sample divided into combinations (a) with and (b) without 
a Df consistent with D* s {2112) + — » Df'y decay. The curves 
and shaded regions, as described in Fig. 7, represent the result 
of a likelihood fit. 



that the curve cannot be distinguished in Fig. 7. 

The -D s j (2460) + reflection arises from contamination 
from D*(2112) + decay It is an interesting exercise to 
identify and separate some of this background. This can 
be accomplished by searching for any 7 candidates that, 
when combined with the Df in the same event, produce 
a L>+7 mass within 15 MeV/c 2 of the D*(2112)+ mass. 
Those Dfir combinations in which such a match is not 
found will contain a smaller proportion of -D s j (2460) + re- 
flection, whereas the remaining Dfir combinations will 
contain fewer D*j(2317) + decays. This is indeed the case 
as illustrated in Fig. 9. The same likelihood fit procedure 
used for the entire sample is repeated for these Dfir sub- 
samples, including the MC prediction of the yield and 
shape of the D s j(2460) + reflection. The fit results are 
consistent with the data. 

As can be seen in Figs. 7 and 9, the £>* J (2317)+ 
line shape derived from MC simulation and used un- 
changed in the likelihood describes the data well. Since 
the MC simulation is configured with an intrinsic width 
(0.1 McV) nearly indistinguishable from zero, it fol- 
lows that the data are consistent with a zero width 
D* sJ (2317)+ meson. 

To extract a 95% CL upper limit on the intrinsic 
£)* J (2317)+ width, the D* sJ (2317)+ line shape is con- 
volved with a rclativistic Breit-Wigncr function aj(m) 
with constant width T: 



aj(m) 



m T 



(to 2 , — m 2 ) + TOqT 2 



(11) 



The fit is then repeated in its entirety at incremental 
steps in T to produce a likelihood curve. Integrating this 
curve as a function of T produces a 95% CL upper limit 
of r < 1.9 MeV (statistical error only). 

In order to produce an estimate of D*j(2317) + yield, 
the fit results must be corrected for selection efficiency. 
This efficiency is calculated using a D*j(2317) + MC sam- 
ple and is p* dependent. Since the p* distribution ob- 
served in data does not exactly match the MC simulation, 
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TABLE III: A summary of systematic uncertainties for the 
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FIG. 10: Corrected D* 7 (2317)+ -» D+tt yield as a function 
ofp*. 



it is important to take into account this p* dependence. 
Two methods are used to do this. The first is to weight 
each D+tt° combination by the inverse of the selection 
efficiency before applying a likelihood fit. After correct- 



ing for absolute D+ 



N(D* sJ (2317) + - D" 3 



yield (Section VI) , the result 



,Z?+^0tt+) = 26 290 ±650 



for p* > 3.2 GeV/c (statistical error only). 

The second method is to divide the D+tt° sample into 
bins of p* . A likelihood fit is applied to each bin and the 
yield corrected for the average selection efficiency in that 
bin. The result is the p* distribution shown in Fig. 10. 
The total yield from this method is 26470 ± 660 (statis- 
tical error only). 

The systematic uncertainties for the D* sJ (2317) + mass 
and yield are summarized in Table III. The uncertainties 
in D s ,/(2460) + yield are calculated in the same fashion. 
The assumed D~f mass value and the 1% relative un- 
certainty in the EMC energy scale are the two largest 
contributors to the error on the D*j(2317) + mass. Un- 
certainties in the signal shape produce the largest uncer- 
tainties in £>* 7 (2317) + yield and width. 

For example, the amount of D s j(2460) + reflection is 
proportional to the D s j(2460) + yield, which, as will 
be discussed in Section XI, has an 11% uncertainty. 
Adjusting the contribution of this reflection in the fit 
with no other changes produces a relative uncertainty 
of 2.0% in the D* 7 (2317)+ yield with little change in 
the L>* 7 (2317)+ mass and width limit. If the likeli- 
hood fit is allowed to choose a Z? s ,/(2460) + reflection 
yield that best matches the data, little change in either 
the .D*j(2317) + mass or yield is observed. The limit 
on the intrinsic width, however, increases to 3.4 MeV, 
since reducing the -D s j(2460) + reflection allows the ob- 
served mass line shape to accommodate a larger intrinsic 
£I* J (2317) + width. 



# s *j(2317) H 
final state. 



mass and yield from the analysis of the D a 





Mass 


Relative 


Source 


(MeV/c 2 ) 


yield (%) 


Df mass 


0.6 


— 


EMC energy scale 


1.3 


— 


D sJ (2460) + reflection size 


<0.1 


2.0 


D s ,/( 2460) + mass 


0.1 


0.7 


Detector resolution 


<0.1 


3.2 


D*(2112)+ reflection model 


<0.1 


0.4 


D~f efficiency 


— 


3.9 


7T° efficiency 


— 


3.0 


p* distribution 


— 


0.6 


Df -> <j>ir + yield 


— 


3.9 


Quadrature sum 


1.4 


7.4 



Another uncertainty that has a similar effect is the as- 
sumed Dfn° mass resolution. Small variations in resolu- 
tion, consistent with comparisons of data and MC simula- 
tion of other known particles, can change the D*j(2317) + 
line shape sufficiently to lower or raise yields by 3.2%. 
Allowing better reconstructed resolution provides more 
room for a large intrinsic width, raising the 95% CL for 
T to 3.0 MeV. 

Other uncertainties in the D*j(2317) + yield include 
the accuracy (±3%) of the MC prediction of ir° efficiency, 
the difference of the two methods for correcting for p*- 
dependent efficiency, and the Df — > </)7r + branching frac- 
tion. The total systematic uncertainty is calculated from 
the quadrature sum of all sources. The result is the fol- 
lowing D* J (2317)+ mass: 

to = (2319.6 ± 0.2 (stat.) ± 1.4 (syst.)) MeV/c 2 , 



and the following yields: 



^(25^(2317)- 
a(D 8j {2460y 



D+n°) 
D+ir°) 



(115.8 ±2.9 ±8.7) fb 
(-1.0 ±1.4 ±0.1) fb, 



where the first error is statistical and the second is sys- 
tematic. 

As can be seen in Table III, the determination of 
the D*j(2317) + mass is limited by the understanding 
of the EMC energy scale. A more primitive calculation 
of this energy scale was used in a previous estimate of 
the .D*j(2317) + mass from this collaboration [4], result- 
ing in an mass estimate that is 2.3 MeV/c 2 lighter then 
the estimate presented here. The associated systematic 
uncertainty in this previous work was also incorrectly cal- 
culated. The central value and systematic uncertainty in 
mass reported here reflects the current best understand- 
ing of these calibration issues. 

For the sake of simplicity, in order to incorporate sys- 
tematic effects into a limit on T, the least strict limit 
obtained from the various systematic checks is quoted. 
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FIG. 11: Likelihood fit results for a D*j(2317) + meson of 
instrinsic width (dashed line) r = and (solid line) F = 
3.8 MeV. Shown in comparison is the mass distribution (solid 
points) for the Dfir candidates. 



This produces a 95% CL of T < 3.8 McV. The lineshape 
produced by this limit is illustrated in Fig. 11. 



X. THE D+7 FINAL STATE 

The -D + 7 mass distribution using 7 candidates with 
loose (150 MeV) and final (500 MeV) minimum energy 
requirements is shown in Fig. 12. Some structure in the 
vicinity of the D* 7 (2317) + and D sJ (2460) + masses be- 
comes apparent once the tighter energy requirements are 
applied. The looser requirement is useful for studying 
the Z?*(2112) + peak. A fit to that peak consisting of two 
Gaussians on top of a polynomial background function 
results in a peak D*(2112)+ mass of 2113.8 MeV/c 2 and 
a yield of 75 000 decays, both with negligible statistical 
uncertainties. 

An unbinned likelihood fit is applied to the final -D + 7 
mass distribution in order to extract the parameters 
and yield of the D s j (2460) + meson and upper limits on 
D*j(2317) + decay. For simplicity, the fit is performed 
only for masses between 2.15 and 2.85 GeV/c 2 . The like- 
lihood fit includes five distinct sources of -D + 7 combina- 
tions: 



• D sJ (2460)+ 
. D* sJ (2317)+ 



D + 7 decay. 

D + 7 decay (hypothetical). 



• A reflection from L>* J (2317)+ -*■ D+tt decay in 
which only one of the 7 particles from ir° decay is 
included. 

• A reflection from D Si7 (2460)+ -* D*(2112) + tt° de- 
cay in which only one of the 7 particles from tt° 
decay is included. 

• Background from both unassociated Df and 7 
mesons and the high-mass tail from D*(2112) + — ► 
D+7 decay. 
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FIG. 12: The invariant Df'y mass distribution using 7 candi- 
dates with loose (150 MeV) and final (500 MeV) energy re- 
quirements. The insert focuses on the region of the Dg (2112) + 
meson using the loose energy requirement. The curve repre- 
sents the fit described in the text. 



The PDF used to describe the mass distribution of each 
of these sources is described below. 

As shown in Fig. 13(b), the reconstructed mass distri- 
bution of the D s j (2460)+ — > Dfj decay, as predicted by 
MC, has a long, low mass tail. To describe this shape, 
the MC sample is fit to a modified Lorentzian function 
F L2 (m): 



F L2 (m) 

S 



«'S 



|l + a A 5 + CL58 2 + 
(l + (5 2 ) a7 
(to - 01) /02 • 



(12) 



The fit results are shown in Fig. 13(b). A similar pro- 
cedure is used for the hypothetical Z?*j(2317) + — ► -D + 7 
decay (Fig. 13(a)). 

Ignoring resolution effects, the £>*j(2317) + reflection 
produces an invariant Z? + 7 mass distribution up to 
a maximum of approximately 140 MeV/c 2 below the 
D*j(2317) + mass. Candidate selection requirements pro- 
duce a distribution that peaks at this limit. Resolution 
effects smear this sharp peak producing the shape shown 
in Fig. 13(c), as predicted by MC. This distribution can 
be reasonably described (in the mass range of interest) 
by a bounded quadratic function convoluted with a dou- 
ble Gaussian. The D s j(2460) + reflection has a similar 
behavior near the D* J (2317)+ mass (Fig. 13(d)). Both 



,(2317)" 



D+7 de- 



distributions overlap the direct D* sJ y 
cay. 

The following function D(m) is used to represent the 
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FIG. 13: The reconstructed Dj~7 invariant mass spectrum 
from MC samples for (a) L>*j(2317)+ and (b) D sJ (2460) + 
decay and (c) D*j(2317)+ -> D+n° and (d) D s ./(2460) + -> 
DJ(2112) + 7r° reflections. The curves are the fit functions de- 
scribed in the text. The D;!j(2317) + signal shapes from (a) 
and (b) are shown for comparison in (c) and (d) in gray. 
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FIG. 14: An example likelihood fit to the D^~7 invariant mass 
distribution. The solid points in the top plot are the mass 
distribution. The open points are the Df sidebands, scaled 
appropriately. The bottom plot shows the same data after 
subtracting the background curve from the fit. Various con- 
tributions to the likelihood fit are also shown. 



remainder of the Df"/ distribution: 

D(m) = l + a 1 cxp[-(m-a 3 ) 2 /a 2 ] ■ (13) 

This includes combinatorial background along with any 
tail from D*(2112)+ -*■ D+-f decay. The MC simulation 
fails to reproduce the shape of this background, cither 
due to unknown cs contributions (for example, higher 
resonant states) or unexpected behavior of the tail of the 
distribution from Z?*(2112) + decay. This issue, combined 
with the complex shapes associated with the D*j(2317) + 
and D s ,/(2460) + reflections, leads to considerable system- 
atic uncertainty in the fit. Likelihood fits under several 
different conditions are attempted in order to understand 
the uncertainty. 

One fit that produces a good representation of the data 
is shown in Fig. 14. In this fit, all parameters except the 
upper mass limit of the D*j(2317) + reflection are allowed 
to vary. The estimated raw £> sJ (2460)+ (D* 7 (2317)+) 
yield from this fit is 920 ± 60 (-130 ± 130). The fitted 
£> sJ (2460)+ mass is (2459.5 ± 1.2) McV/c 2 (statistical 
errors only). 

The D s j(2460) + signal is far enough removed from 
the D* 7 (2317)+ and L> sJ (2460) + reflections that accu- 
rate mass and yield results are obtained. The same two 
methods described in the previous section are used to 
estimate the £> s j(2460) + yield. The first method, us- 
ing p* dependent weights proportional to the inverse of 
efficiency, produces a corrected yield of: 
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FIG. 15: Corrected L> sJ (2460)+ -+ D+7 yield as a function 
of p* . 



for p* > 3.2 GeV/c (statistical error only). The second 
method produces the p* spectrum shown in Fig. 15. The 
total yield from the spectrum is 3080 ± 240 (statistical 
error only), approximately 6.2% lower than the first es- 
timate. 

Systematic uncertainties associated with the assumed 
PDFs for the background and reflection shapes are ex- 
plored using different variations of the likelihood fit. 
Among the variations applied is an alternate description 
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TABLE IV: A summary of systematic uncertainties for the 
Z) s j(2460) + mass and yield from the analysis of the Dfy 
final state. 





Mass 


Relative 


Source 


(MeV/c 2 ) 


yield (%) 


D^t mass 


0.6 


— 


EMC energy scale 


3.7 


— 


D*. 7 (2317) + mass 


0.1 


0.1 


D*j(2317) + reflection shape 


0.1 


1.9 


Detector resolution 


<0.1 


3.5 


Background shape 


0.5 


3.3 


Df efficiency 


— 


3.9 


7 efficiency 


<0.1 


1.8 


p* distribution 


— 


6.3 


Df — > (t>ir + branching fraction 


— 


3.9 


Quadrature sum 


3.7 
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FIG. 16: (a) The sample of Dfir '7 candidates shown in solid 
points. The fit described in the text is overlaid. The open 
points are those candidates which fall in a restricted Df^y 
mass range, (b) The reconstructed Dfn "/ invariant mass 
spectrum from a D sJ (2460) + -> D*(2112) + 7r° MC sample. 
The fit function of Eq. 15 is overlaid. 



of the background shape: 

D'(m) = l + aiexp[-5(TO-a 3 )] + a 2 /(TO - a 3 ) . (14) 

In addition, the MC predictions for the size and shape 
of the D*j (231 7) + reflection are used unaltered (despite 
producing a fit of inferior quality). Large variations of 
raw D*j(2317) + yield of up to 490 events are observed. 

The shape of the D s j(2460) + — > Dfj signal is sensi- 
tive to several factors that are difficult to simulate ex- 
actly, including EMC energy resolution. Variations in 
the assumed resolution are used to study the associated 
systematic uncertainty in yield and mass. The result is 
an uncertainty of 3.5% in yield and no significant change 
in mass. 

All systematic uncertainties for the D s j(2460) + mass 
and yield are listed in Table IV. The result is the follow- 
ing Z) s j(2460) + mass: 

m = (2459.5 ± 1.2 (stat.) ± 3.7 (syst.)) MeV/c 2 , 



and the following yields: 



a(D* sJ (2317)+ 
a(D sJ {2460)+ 



D+j) = (-2.4±2.3±8.9)fb 

L»+ 7 ) = ( 14.4 ±1.0 ±1.4) fb, 



where the first error is statistical and the second is sys- 
tematic. 



XI. THE £>+tt°7 FINAL STATE 



The invariant mass spectrum for all selected Dfir ^ 
candidates is shown in Fig. 16(a). A D s j(2460) + signal 
is apparent. The shape of this signal is characterized by 
applying the following modified Lorcntzian fit function 



F L3 (m) 
S 



= a 3 



1 + <24<5 + 05 tan 1 5(5 



(1 + S 2 ) a6 
(to - 01) /02 • 



(15) 



to a £> sJ (2460)+ -► D* s (2112)+tt° MC sample 
(Fig. 16(b)). A binned x 2 fit to the Dfir°j spectrum 
that includes this shape along with a polynomial de- 
scription of the background produces a D s j(2460) + mass 
(2459.5±2.0) MeV/c 2 and a yield of 560±80 events (sta- 
tistical errors only). 

In the following, it is assumed that the £> s j(2460) + 
meson decays to Dfir 0r y entirely through either of the 
two kincmatically allowed sub-resonant decay modes: 



£> sJ (2460)- 



D*(2112) + Ti a 
D* sJ (2317)+ 1 



D. 



-7T 7 



(16) 



Due to a kinematic accident, the phase space of these two 
sub-resonant modes overlap, as illustrated in Fig. 17. It 
is therefore possible to remove background while retain- 
ing both sub-resonant decay modes by selecting £>+7r°7 
candidates in cither a restricted range of Dfir mass or 
a restricted range of Df-f mass. This analysis uses a re- 
quirement that the Dfj mass must be within 20 MeV/c 2 
of the £>*(2112)+ mass. As shown in Fig. 16(a), the re- 
sulting _D+7r°7 mass distribution in this Dfj signal win- 
dow is considerably cleaner. 

The -D+7 requirement introduces a source of back- 
ground that peaks underneath the D s j(2460) + signal. 
This background is a reflection from D* sJ (231 7) + — ► 
-D+7T decays that are not associated with any 
D s j(2460) + decay. The reflection arises because, as il- 
lustrated in Fig. 17, any D*j(2317) + signal that is com- 
bined with a 7 candidate that produces a Dfj mass 
near the D*(2112) + meson results in a Dfir j mass near 
the D s j(2460) + meson. If the Df-f mass requirement is 
shifted upwards or downwards, this -D*j(2317) + reflec- 
tion shifts up and down in Z?+7r°7 mass by a predictable 
amount. 

Another type of background that behaves similarly to 
the D*j(2317)+ reflection is L> s j(2460)+ -► D+tt 7 de- 
cay in which the wrong 7 candidate is chosen. This 
type of background is slightly wider and smaller than the 
£>*j(2317) + reflection, but otherwise has a similar mass 
distribution. To describe this background contribution, it 
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D sj (2460)+ -> D* n° 7 



DL(2317)* 




--- D;(2112) + 
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D^tt Mass (GeV/c 2 ) 



FIG. 17: The light gray region indicates the range of D^n 
and D^7 mass that is kinematically allowed in the decay of 
an object of mass 2458.0 MeV/c 2 to D+7r°7. The lines mark 
the kinematic space associated with decays which proceed 
through an intermediate Dj(2112) + or D*j(2317) + meson. 



is assumed that the D s j(2460) + decays entirely through 
£>*(2112) + 7T and is produced at a rate comparable to 
previous measurements. Both assumptions need not be 
entirely accurate since this background has a relatively 
small contribution. 

and £> sJ (2460)+ re- 
mass selection win- 
2 away from 



To characterize the D* 7 (2317) 



flections, upper and lower D+7 
dows are chosen centered at ±60 MeV/c 
the D*(2112) + mass. The mass distribution from MC 
samples of L>* 7 (2317) + -> D+tt° and D sJ (2460)+ -► 
D*(2112) + 7 decay are shown in Fig. 18 for the signal 
and two sideband Dfj mass windows. The shape of the 
two combined reflections in all three cases can be suc- 
cessfully described by a fit to a Gaussian. 

To determine the mass, width, and yield of the 
.D s j(2460) + meson, an unbinncd likelihood fit is applied 
to the D+7r°7 mass distribution of candidates selected in 
the Df"f signal window. This fit includes the following 
contributions: 

• D sJ (2460)+ -> £>+tt°7 decay. 

• The combined reflections from D* 7 (2317) + — ► 
L>+7T° decay and from D Si7 (2460)+ -> L»*(2112)+7 
decay in which the incorrect 7 candidate is chosen. 



A reflection from D*(2112)" 



D+7 decay in 



which an unassociated 7 candidate is added to form 
a fake tt° candidate. 

• Smooth background sources that do not have any 
peaking behavior. 

The D*(2112)+ -> D+7 reflection is similar to that ob- 
served in Dfir combinations (see, for example, Fig. 7). 
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FIG. 18: The combined Dfiv j invariant mass distribution 
(solid points) as obtained from MC samples for the com- 
bination of the D*j(2317)+ reflection and D sJ (2460)+ -+ 
Da(2112) + 7r° decays in which the incorrect 7 is chosen. 
The D s j (2460) + contribution alone is shown in open points. 
Shown are the (a) upper, (b) signal, and (c) lower _D7J~7 mass 
selection windows. The curves are fits to Gaussian distribu- 
tions. 



The smooth background is represented by the C(m) func- 
tion described in Eq. 6. 

Two similar fits excluding the D s j(2460) + signal are 
applied to the upper and lower Dfj mass samples. These 
fits suggest that the MC prediction of the absolute rate of 
the D*(2112) + — > Dfj reflection is approximately 21% 
too low. The fit models are adjusted accordingly. The fit 
results for all three -D+7 mass ranges after this adjust- 
ment are shown in Fig. 19. The result is a D s j(2460) + 
mass of (2458. 6± 1.3) MeV/c 2 and a raw yield of 560±40 
events (statistical errors only). 

Although the overall size of the D* ,-(2317)+ -» D+n" 
reflection is allowed to vary, the MC prediction for the 
relative contributions in the three Dfj mass windows is 
preserved. Since the size of this reflection is adequately 
modeled in the two Dfj sidebands, there is some confi- 
dence that the size is well established in the D+7 signal 
window. 

Note that the size of the smooth background is rel- 
atively larger in the signal Dfj window due to contri- 
butions from D*(2112) + — > Dfj decay. In addition, 
if there were significant non-resonant contributions to 
D s j(2460) + decay, peaks at the D s j(2460) + mass would 
be visible in the two Dfj sidebands. No such evidence 
is visible. 

The two methods described in the two previous sec- 
tions are used to estimate the D s j(2460) + yield. The 
first method, using p* dependent weights proportional to 
the inverse of efficiency, produces a corrected yield of: 



N(D sJ (2A60) + ^Dtw 7 ,D^ 



</>7T H 



9690 ± 790 



for p* > 3.2 GeV/c (statistical error only). The second 
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D^ 7T 7 Invariant Mass (GeV/c ) 

FIG. 19: The invariant mass distribution of D^tv°j candi- 
dates in the (a) upper, (b) signal, and (c) lower D^"7 mass 
selection windows for (solid points) the D+ signal and (open 
points) Df sideband samples. The curves represent the fits 
described in the text. The dark gray (light gray) region cor- 
responds to the predicted contribution from the D*j(2317) + 
(£>*(2112)+) reflection. 



method produces the p* spectrum shown in Fig. 20. The 
total yield from the spectrum is 9890 ± 810 (statistical 
error only). 

The systematic uncertainties in the mass and yield of 
the jD s j(2460) + meson are shown in Table V. As de- 
scribed previously, the size of the D*j(2317) + reflection 
was adjusted in the fit to match the Dfj sideband sam- 
ples. If the size of the reflection is taken unchanged from 
MC predictions, the D s ,/(2460) + yield increases by 7.3%. 
A second likelihood fit described later in this section used 
to distinguish between the two D s j(2460) + sub-resonant 
decay modes also produces an estimate of D s j(2460) + 
yield. The difference between the two fits is treated as 
a systematic uncertainty. The total D s ,/(2460) + yield, 
without distinguishing between the two possible sub- 
resonant decay modes, is measured to be 

ct(L» sJ (2460) + -► £>+7r°7) = (42.7 ± 3.5 ± 4.2) fb, 

where the first error is statistical and the second is sys- 
tematic. The complete D s j (2460) + mass result is: 

m = (2458.6 ± 1.0 (stat.) ± 2.5 (syst.)) MeV/c 2 . 

The LW(2460) + signal PDF used in the likelihood 
fit of Fig. 19 includes a D s j(2460) + intrinsic width of 
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FIG. 20: Corrected L> sJ (2460)+ -> D+it ^ yield as a function 
of p* . 



TABLE V: A summary of systematic uncertainties for the 
D s j(2460) + mass and yield from the analysis of the D+7T 7 
final state. 





Mass 


Relative 


Source 


(MeV/c 2 ) 


yield (%) 


D+ mass 


0.6 


— 


EMC energy scale 


2.4 


— 


D*j(2317) + reflection size 


0.3 


7.3 


D* aJ (2317) + mass 


0.1 


1.2 


Detector resolution 


<0.1 


2.1 


D*(2112)+ reflection model 


— 


— 


Fit method 


— 


1.4 


D^ efficiency 


— 


3.9 


7T° and 7 efficiency 


— 


3.0 


p* distribution 


— 


2.1 


Df -> 07T+ yield 


— 


3.9 


Quadrature sum 


2.5 


10.2 



r = 0.1 McV. Larger intrinsic widths do not result 
in any significant improvement of the fit. After apply- 
ing the same likelihood-integration technique described 
in Section IX for the Dfir final state, and including the 
systematic effects listed in Table V, the result is a 95% 
CL limit of T < 6.3 MeV. 

Having established a Z) s ,/(2460) + — ► Dfir "/ signal, it 
is now necessary to distinguish between the two possible 
sub- resonant decay modes shown in Eq. 16. These two 
decay modes can be distinguished by their Dfir and 
.D+7 invariant mass distributions, as shown in Fig. 21. 
The distributions for the D sJ (2460) + -> £>*(2112)+7r° 
subresonant mode are determined using a MC sample. 
The reconstructed Dfj mass distribution, which is rel- 
atively narrow (as it arises from D*(2112) + decay), is 
represented by a \ 2 fit to the F^ function (Eq. 15). The 
wider Dfir mass distribution is accurately modeled by 
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FIG. 21: The reconstructed (a) D+tt and (b) D+7 in- 
variant mass distributions for the two possible D s j (2460) + 
sub-resonant decay modes. The distributions from a 
D sJ (2460) + -> £)*(2112) + 7T° MC sample are shown in points. 
The curves are the fits described in the text. The shaded re- 
gions are the shapes assumed for D s j(2460) + — » D*j(2317) + 7 
decay. 



Both 



a square function smeared by a double Gaussian, 
fits are shown in Fig. 21. 

In contrast, for the D sJ (2460) + -> Z?* J (2317)+7 sub- 
resonant decay mode, the Dfir mass distribution is nar- 
row and the Dfj mass distribution is wide. The D+ir° 
distribution is determined using a D* 7 (2317) + — ► Dfw° 
MC sample. The D+j distribution is calculated using the 
parameters determined from the D+ir° distribution from 
D s j(2460) + -> D*(2112)+7r° described above converted 
to the appropriate kinematic range. The shapes assumed 
for both D sJ (2460)+ -> L>* 7 (2317)+7 mass projections 
are shown in gray in Fig. 21. 

The Dfir and Dfj mass distributions of the sig- 
nal cannot be explored without correctly subtracting 
backgrounds from unassociated D*(2112) + — ► D+7 and 
D*j(2317)+ -> D+tt° decay. This subtraction is per- 
formed by a two-dimensional unbinncd likelihood fit ap- 
plied to the D+ir° and Dfj mass distributions of the 
data. The likelihood fit is restricted to the data sam- 
ple contained inside the grid shown in Fig. 22. This fit 
includes five sources of D+7r°7 candidates: 

• Combinatorial background represented by a two- 
dimensional quadratic function. 

• Background from D*(2112) + — ► Dfj decay com- 
bined with unassociated ir° candidates represented 
by a Z?*(2112) + line shape in the -D+7 mass and as 
a linear function in Dfir mass. 

• Background from D* j (2317) + — ► Dfn° decay com- 
bined with unassociated 7 candidates represented 
by a D*j(2317) + line shape in the Dfn° mass and 
as a linear function in Df"/ mass. 

• A signal from £> sJ (2460) + -» L>*(2112)+7r° with 
Dfir and Dfj mass distributions represented by 



the curves in Fig. 21. 

A signal from D sJ (2460)' 



D* sJ (2317)+j with 



2.3 



2.2 - 



2.1 - 




FIG. 22: The Df'j versus D+ir° mass distributions for 
the D+7r°7 candidates. The horizontal (vertical) band 
corresponds to background from D*(2112) + — > Df'y 
(D*j(23l7) + — ♦ Dfiv ) decay. The excess of candidates near 
the crossing of these two bands is the D s j (2460) + signal. 
The curve indicates the region of phase space in which the 
D s j(2460) + decay is kinematically restricted. The grid iden- 
tifies the subsample of candidates used in the likelihood fit 
shown in Fig. 23. 



The result of this likelihood fit is shown in Fig. 23, 
divided into the regions delineated by the grid shown 
in Fig. 22. The fit produces an adequate model of the 
data in all regions. The result (statistical errors only) is 
a total yield of 520 ± 50 L> s ./(2460)+ -> D+tt -/ decays 
with a fraction of (2.5 ± 8.8)% proceeding through the 



Dfir and Dfj mass distributions represented by 
the gray regions in Fig. 21. 



£>*j(2317) + 7 channel, the former number being some- 
what smaller than the yield determined by the D3 L 7r°7 
mass fit (Fig. 19), though consistent within systematic 
uncertainties. Based on these results, it appears that 
the decay -D s j(2460) + — > D+7r°7 can be successfully 
described as proceeding entirely through the channel 
D*(2112)+7T°. 

The systematic uncertainties listed in Table V can be 
applied to the above fit results. In combination with 
the results of the fit to the Z?+7r°7 mass distribution of 
Fig. 19, and treating correlated systematic uncertainties 
in the appropriate fashion, the following yields for the 
subrcsonant specific decays are obtained: 

a(D sJ (2A60)+ -> D*(2112)+tt°) = (41.6 ± 5.1 ± 5.0) fb 
ct(D sJ (2460)+ -* -D^(2317)+7) = ( 1.1 ± 5.1 ± 5.0) fb, 

where the first error is statistical and the second is sys- 
tematic. 

A simple helicity analysis is used to test the J p as- 
signment of the -D s j(2460) + meson under the assump- 
tion that the decay Z? s ,/(2460) + — ► Z?+7r°7 proceeds en- 
tirely through the subrcsonant mode D*(2112) + 7r°. This 
analysis is performed in terms of the helicity angle #/j, 
defined as the angle of the 7 in the D*(2112) + centcr- 
of-mass frame with respect to the D*(2112) + direction. 
Since the D*(2112) + is a vector particle, the helicity dis- 
tribution must consist of some combination of the zero 
helicity distribution Hq: 



H = sin d h 



(17) 
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FIG. 23: The Dfn and -Dj~7 invariant mass distributions 
for Z?j"7r°7 candidates that fall within the indicated por- 
tions of the grid shown in Fig. 22. The histograms represent 
the results of a likelihood fit. The light gray region corre- 
sponds to combinatorial background. The medium gray (dark 
gray) region represents the fitted fraction of D s ./(2460) + — » 
D* 3 (2112)+ir° (D sJ (2460)+ - L>: J (2317)+ 7 ). 



TABLE VI: Helicity distributions in the decay D sJ (2460)+ -» 
D*(2112) + 7T° for hypothetical D 3 j(2460) + spin-parity assign- 
ments. 



J' 






Helicity Distribution 



i + 



2-,3+,4- 



(decay is forbidden) 

Hi 

any combination of Hq and Hi 



and the helicity = ±1 distribution H\: 



1 



fli = 2 (1 + cos 2 $ h ) 



(18) 



As listed in Table VI, the expected combination of H 
and Hi depends on the assumed D s j(2460) + spin and 
parity. 

To measure the helicity distribution, the Dfir°j can- 
didates are divided into five bins of cosi?/,. The D+7r°7 
mass fit of Fig. 19 is repeated for each of these subsam- 
ples using p*-depcndent weights inversely proportional 
to the selection efficiency in order to correct for accep- 
tance. The result is shown in Fig. 24. The integral of the 
following function is calculated in each cosfih bin: 
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FIG. 24: EfBciency-corrected yields in five cos'&h bins for the 
decay D sJ (2460)+ -> D*(2112) + 7r°. The solid histogram is 
the result of a fit to the function described in the text. The 
dashed histogram is a similar fit with the a-z parameter fixed 
to zero. 



A x 2 fit is used to determine the most likely value of a 2 . 
The result is a 2 = 0.76 ± 0.14 (statistical errors only). 

The same procedure can be repeated after each rele- 
vant systematic check listed in Table V is performed. The 
differences in a 2 values so obtained are added in quadra- 
ture to estimate the total systematic uncertainty The 
final result is: 

a 2 = 0.76 ± 0.14 (stat.) ± 0.06 (syst) , (20) 

where a 2 = (a 2 = 1) corresponds to a D s j(2460) + of 
helicity zero (±1). This value of a 2 deviates from zero 
by 5.1 standard deviations, which strongly disfavors the 
J p = 0~ interpretation of the D s j(2460) + while remain- 
ing consistent with a J = 1 or higher interpretation of 
cither parity 



XII. THE Dfir°ir° FINAL STATE 



The D 



+ 7T 7r 



final state contains potential contribu- 



tions from L> s j(2460) + -» D*(2112)+7r° subresonant de- 
cay through the channel £>*(2112) + — > Z?+7r°, which has 
a branching ratio of 5.8% [15]. Since this sub-resonant 
mode is more efficiently investigated using the Dfir 0r Y 
final state (as discussed in the previous section), it is re- 
moved from the D+7r°7r° sample by requiring the Dfir 
invariant mass to be greater than 2117.4 MeV/c 2 for both 
7T° candidates. This requirement excludes the edges of 
the .D*j(2317) + and -D s j(2460) + phase spaces, as illus- 
trated in Fig. 25. This figure also demonstrates how little 
phase space is available to the D*j(2317) + meson in this 
decay in comparison to the Z) s j (2460) + meson. 

The invariant mass distribution of the selected 



D 



+7r°7r° 



F h =a 1 [(l-a 2 )H Q +a 2 Hi} 



(19) 



candidates is shown in Fig. 26. There is no ev- 
idence of £>*j(2317) + or D s ,/(2460) + meson decay, nor 
is there evidence of any structure in the background. 
The mass distribution is fit using a likelihood function 
that consists of the smooth background function C(m) 
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FIG. 25: The Dalitz phase space available to (a) the 
D*j(2317) + and (b) the D sJ (2460) + mesons in the D+ttV 
final state. The dashed horizontal line corresponds to the 
D*(2112) + mass. The dark shaded regions are those parts 
of the phase space removed by the requirement m(Z)+ 7r°) > 
2117.4 MeV/c 2 . The dotted curve in (b) is the D*j(2317) + 
phase space drawn for comparison. 
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FIG. 26: The invariant mass distribution for (solid points) 
D^7r°7r° candidates and (open points) the equivalent using 
the Df sidebands. The curve represents the likelihood fit 
described in the text. The insert focuses on the low mass 
region. The dotted line in the insert indicates the D*j(2317) + 
mass. 



(Eq. 6) and D* 7 (2317) + and D sJ (2460)+ meson contri- 
butions, the latter two having yields that are allowed to 
fluctuate to negative values. The shape of the signals is 
modeled by double Gaussians, the parameters of which 
are determined by fits to MC samples. The masses of the 
D* 7 (2317)+ and L> sJ (2460)+ mesons are fixed to values 
measured in the previous sections. 

The result of the likelihood fit is shown in Fig. 26 
and produces a D* 7 (2317)+ (L> s ./(2460)+) raw yield of 
0.2 ±3.9 (5 ±10) (statistical errors only). Efficiency cor- 
rections are calculated by applying a Z?+7r°7r° candidate 
weight that is inversely proportional to p*-dependent se- 
lection efficiency, calculated from MC samples. This pro- 
cedure produces corrected D*j(2317) + and D s j(2460) + 
meson yields of 1900± 1600 and 1300± 1200, respectively 



(statistical errors only). 

Various systematic uncertainties are considered. Ac- 
cording to MC simulation, the selection efficiency varies 
by as much as 20% across the Dalitz plot. Since no 
specific Dalitz distribution is assumed, this variation is 
translated directly into a multiplicative systematic un- 
certainty. An alternate shape C^(m) for the background 
is considered: 



C' 2 (m) 



-«r(l 



a2rn + a^m 



(21) 



A fit with this background increases the raw D*j(2317) + 
yield by 0.5 events. Detector resolution and D*j(2317) + 
and .D s j (2460) + meson mass variations have similar ef- 
fects. The final results are 



a(D* sJ (2317)+ 
a(£> s j(2460)+ 



L>+7T 7r°) = (8.7±6.9±5.0)fb 
L>+7T 7T°) = (5.5 ±5.4 ±2.4) fb, 



where the first error is statistical and the second is sys- 
tematic. 



XIII. THE £>+77 FINAL STATE 

It is assumed that there are three possible contribu- 
tions to the decay of the D* 7 (2317)+ and £> sJ (2460) + 
mesons to the -D+77 final state: 



Two body decay to Dfir followed by 7T° 



77. 



• Subrcsonant 
£>*(2112)+- 



decay 

Dfl- 



to D*(2112)+7 followed by 



• Non-resonant decay directly to .D+77. 

The -D+7T final state, already studied in Section IX, is 
removed from the Dfjj sample by the 7 selection re- 
quirements. Potential background from ir° mass tails is 
removed by further requiring the 77 invariant mass to be 
less than 100 MeV/c 2 or greater than 170 MeV/c 2 . The 
remaining two Dfjj contributions are treated by divid- 
ing the sample into two portions, one rich in Z?*(2112) + 
decay (referred to as the D*(2112) + sample), and the 
remainder (referred to as the non-resonant sample). 
A Dfjj candidate is placed in the £>*(2112) + sam- 
ple if cither 7 produces a Dfj invariant mass within 
15 MeV/c 2 of the PDG value for the D*(2112)+ mass 
(2112.4 MeV/c 2 ) [15]. 

Figure 27 illustrates how the the Df-f-f candidates are 
divided in terms of the phase space of D*j(2317) + and 
D s j(2460) + meson decay. The D+77 invariant mass dis- 
tribution of the two samples is shown in Fig. 28. No clear 
Z?* J (2317)+ or D sJ (2460)+ signal is observed. The mass 
distributions contain structure associated both with the 
7T° veto requirements and with the D*(2112) + sample 
selection requirements. Much of this structure can be 
avoided by applying likelihood fits in a restricted mass 
range. Two types of background that produce structure 
that cannot be so avoided arc described below. 
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FIG. 27: The Dalitz phase space corresponding to the DJ'y-y 
samples for (a) D*j(2317)+ and (b) D sJ (2460)+ mesons. 
The dark shaded and light shaded regions correspond to the 
Dg(2112) + and non-resonant samples, respectively. The hor- 
izontal (vertical) dashed line indicates the D*(2112) + (7r ) 
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FIG. 29: (a) Shown on top as a function of Dfy and Df-y-y 
masses are the bands corresponding to the D*(2112) + sample 
requirement (-D+7 mass within 15 MeV/c 2 of the D*(2112) + 
mass) and the effect of this requirement on the other 7 candi- 
date. The overlap of the two bands produces the shape shown 
below, (b) The D+77 mass distribution of the D sJ (2460)+ -> 
Ds(2112) + 7r reflection from a MC sample. The curve is the 
fit described in the text. 
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FIG. 28: The Dfyy invariant mass distribution (solid points) 
for candidates in (a) the Dg(2112) + sample, (b) the non- 
resonant sample, and (c) either sample. The open points are 
the corresponding distributions for the Df sideband. The ver- 
tical dashed lines indicate the D*j(2317) + and D s j(2460) + 
meson masses. 



Either of the two 7 candidates in each -D+77 combi- 
nation provides two opportunities for the D+77 combi- 
nation to be placed into the D*(2112) + sample. At a 
specific value of D+77 invariant mass, however, if one 
7 candidate falls inside the D*(2112) + mass window, 
the other 7 is likely to do the same, due to kinemat- 
ics. This produces a deficit of candidates at this mass in 
the D*(2112) + sample along with a corresponding excess 
of candidates in the non-resonant sample. The -D^~77 in- 
variant mass distribution of this excess can be approxi- 
mated by the triangular shape centered at 2255.9 MeV/c 2 
shown in Fig. 29(a). 



which only ap- 
1 reflection from 
L> s j(2460)+ -► £>*(2112)+7r° decay in which one 7 par- 
ticle from the 7T° decay is ignored. This reflection is sup- 
pressed but not eliminated by the requirement that each 
7 candidate not belong to a fiducial 7r° candidate. The 
MC prediction of the shape of this reflection is shown 
in Fig. 29(b). The shape is accurately modeled by a 
quadratic polynomial bounded on two sides and smeared 
by a double Gaussian. 

The unbinned maximum-likelihood fit to the D+77 
mass distribution of the D*(2112) + sample is performed 
between masses of 2.22 and 2.70 GeV/c 2 . This mass dis- 
tribution has a cusp near a mass of 2.4 GeV/c 2 that is dif- 
ficult to describe using a simple polynomial. Instead, the 
combinatorial background is parameterized using an em- 
pirical function K{m) composed of a line and a parabola 
that intersect at one point. To make the function smooth, 
a cubic spline Ks{m) is used near the intersection point 
(to = a±) to extrapolate between these two polynomials: 



(12 + a^m m < a\ — 5 

K{m)={ K s {m) \m-ai\<6 

&2 + &4TO, + a^m 2 m > ai + 5 



(22) 



where the value S = 20 MeV/c 2 is chosen to approxi- 
mately match the resolution. 

The reconstructed mass distributions for the hypothet- 
ical decay of D*j (2317)+ and £> sJ (2460)+ mesons to 
Z5+77 and D*(2112) + 7 are modeled by the functional 
form of Eq 12, the parameters of which are determined 
by fits to the corresponding MC samples. 

The unbinned likelihood fit to the -D+77 mass distribu- 
tion of the D*(2112) + sample includes five components: 

1. £>*,-(2317) + -► D*(2112)+7 decay (hypothetical). 

2. £> sJ (2460)+ -► D*(2112)+7 decay (hypothetical). 

3. A deficit of events of the shape described in 
Fig. 29(a) but of variable size. 
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FIG. 30: The .D+77 invariant mass distribution for the 
Z?a(2112) + sample is shown on top. The thick curve is the 
result of the likelihood fit. The mass distribution after sub- 
tracting the portion assigned by the fit to combinatorial back- 
ground is shown on bottom. 



FIG. 31: The Dj~77 invariant mass distribution for the non- 
resonant sample is shown on top. The thick curve is the result 
of the likelihood fit. The mass distribution after subtracting 
the portion assigned by the fit to combinatorial background 
is shown on bottom. 



4. A contribution from the D sJ (2460)+ — ► 
D*(2112) + 7T° reflection, of fixed shape based 
on the fit of Fig. 29(b) and with a yield consistent 
with the results of Section XL 

5. A background described by the function of Eq. 22. 

The result of applying the fit function on the D*(2112) + 
sample is shown in Fig. 30. A raw D*j(2317) + 
(L> sJ (2460) + ) yield of -40 ± 110 (-50 ± 140) is obtained 
(statistical errors only). 

The proximity of the £> sJ (2460)+ -► D*(2112)+7T° re- 
flection to potential D* 7 (2317) + and D sJ (2460) + signals 
coupled with the unknown shape of the combinatorial 
background leads to large uncertainties in the fit results. 
For example, if a third-order polynomial is used in place 
of the function K(m) to describe the smooth background, 
the D* 7 (2317)+ and D sJ (2460)+ raw yields increase by 
260 and 150 candidates, respectively. Although the result 
of this fit is a less faithful representation of the mass dis- 
tribution of the data, as a conservative estimate the entire 
difference is quoted as a systematic uncertainty. Other 
systematic checks performed include a variation of the 
range in D+77 mass over which the fit is applied, consid- 
eration of uncertainties in detector resolution, and varia- 
tions of the D* J (2317)+ and L> sJ (2460) + masses and the 



L>;(2112)+7r° rcficc- 



relative size of the L> sJ (2460) 4 
tion. 

The raw D* J (2317)+ and D sJ (2460)+ yields are cor- 
rected for selection efficiency by weighting each -D+77 
combination by the inverse of the p*-dependent selection 



efficiency, determined using a MC sample, 
for the D*(2112)+7 final state are: 



The results 



a(D* 8j (2317Y 
a(D sJ (2460)- 



D*(2112)+7) 
£> s *(2112)+ 7 ) 



(-0.5 ±3.2 ±8.1) fb 
(-0.9±3.5±4.1)fb, 



where the first error is statistical and the second is sys- 
tematic. 

The D^77 mass distribution of the phase-space sub- 
sample has a cusp that is more pronounced than in 
the .D*(2112) + 7 sample but at a mass well below the 
D*j(2317) + mass. This cusp is avoided by restricting 
the fit to masses above 2.2 GeV/c 2 . In this mass range a 
simple polynomial is sufficient to describe combinatorial 
background. The unbinned likelihood fit includes four 
components: 



1. ££,(2317)- 

2. 7J sJ (2460)- 



-D+77 decay (hypothetical). 
-D+77 decay (hypothetical). 



3. An excess of events of the shape described in 
Fig. 29(a) but of variable size. 

4. Combinatorial background represented by a third- 
order polynomial. 

The result of applying the fit function on the non- 
resonant sample is shown in Fig. 31. A raw D* sJ (231 7) + 
(L> sJ (2460)+) yield of 190 ± 120 (80 ± 160) is obtained 
(statistical errors only). 
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As in the D*(2112)+ sample, raw D* sJ (23l7) + and 
D s j(2460) + yields are corrected for selection efficiency 
by weighting each Dfn combination by the inverse 
of the p*-dcpendcnt selection efficiency. This efficiency 
was determined using a MC sample that simulated the 
non-resonant decays of the D* J (2317)+ and £> s j(2460)+ 
mesons to Dfn using phase space (such that the Dalitz 
plot was evenly populated). To simulate an alternate 
decay model, the MC samples were weighted to form a 
cos 2 $d 7 distribution, where $£> 7 is the angle between 
the D+ and each 7 in the D* sJ (2317) + and L> sJ (2460)+ 
center-of-mass frame. This weighting had the effect of 
reducing the selection efficiency by 25% and 15% for the 
D* sJ (2317) + and L> s ./(2460)+ mesons, respectively. 

Other systematic checks include variations in detec- 
tor resolution, the D*j(2317)+ and D sJ (2460)+ meson 
masses, and the fit range. The final results for the non- 
resonant D+77 final state are: 



<7(D* 7 (2317)^ 
a(D sJ (2460r 



Dfll) = (7.4 ±4.5 ±2.2) fb 
Dill) = (3.5 ± 4.3 ± 1.7) fb, 



where the first error is statistical and the second is sys- 
tematic. 



XIV. THE £»+tt ± FINAL STATES 



The invariant mass distributions of the D+ir~ and 
Dfir + candidates are shown in Fig. 32 for masses be- 
tween 2.25 and 2.61 GeV/c 2 . No structure is apparent 
and both distributions can be adequately described by 
a second-order polynomial. Since no signal is appar- 
ent, the remaining task is to place limits on the yield 
of hypothetical doubly-charged (D*j(2317) ++ ) or neu- 
tral (D* sJ (2317)°) partners of the D* J (2317)+ meson. To 
do so requires a few assumptions: 
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FIG. 32: The invariant mass distributions (solids points) of 
(a) Dfn~ and (b) Dfjr + candidates and (open points) the 
equivalent for the D* sidebands. The shaded regions indicate 
the range of assumed D*j(2317) ++ and D*j(2317)° masses. 
The vertical dotted line marks the D s j(2460) + mass. The 
curves are the fits described in the text that produce the 
largest yield within the shaded region. 



small systematic uncertainties. The largest uncertainty 
arises from the assumed shape of the background and is 
estimated by substituting a third-order polynomial for 
the second-order one. There is also a 1.3% relative un- 
certainty in the reconstruction efficiency of each it can- 
didate. The results from the fits that produce the largest 
yields are: 



The intrinsic width T of either the £>* J (2317)° or a(D* sJ (2317)° -> D 



D* sJ (23\7) ++ particle is too small to be resolved by 
the detector. 

• The mass of the DJ 7 (2317)° or D* sJ {23\7) ++ parti- 
cle is within ±10 MeV/c 2 of the D* sJ (2317) + mass. 

A series of unbinned likelihood fits based on these as- 

and Df-K + mass 
or 



a(D;j(2M7)- 



DtTT+) = 



= (1.07 ± 0.44 ± 0.10) fb 
= (0.74 ± 0.37 ± 0.07) fb, 



sumptions are applied to the Dj 



distributions. Included in these fits is a D* sJ (2317)° 
D*j(2317) ++ signal modeled using a line shape extracted 
from a fit to a D; J (2317)° -» D+tt- MC sample. Ac- 
cording to this fit, the mass resolution is approximately 
1.3 MeV/c 2 . To avoid potential statistical biases, the 
mass of each hypothetical particle is fixed at a specific 
value for each fit. Several such fits are applied with the 
assumed mass placed between 2307.3 and 2327.3 MeV/c 2 
at intervals of 1 MeV/c 2 . The fits that produce the 
largest yields are shown in Fig. 32. 

The two final states discussed in this section involve 
only charged particles and are thus subject to relatively 



where the first error is statistical and the second is sys- 
tematic. These results will be used to calculate upper 
limits on the cross section. 



XV. THE D+n+n- FINAL STATE 

The invariant mass distribution of the D+7r + 7r~ candi- 
dates is shown in Fig. 33. Clear peaks from D s j(2460) + 
and .D s i(2536) + decay are apparent. The distribution 
has no additional structure of any statistical significance. 
To determine the mass and yield of the £> s j(2460) + me- 
son, and to place limits on _D*j (2317) + decay, an un- 
binned likelihood fit is applied to the D+7r + 7r _ mass dis- 
tribution. This fit includes the following contributions: 

• D* J (2317) + -► D+tt+tt- decay (hypothetical). 
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FIG. 33: The invariant mass distribution of (solid points) 
D^7r + 7r~ candidates and (open points) the equivalent using 
the Df sidebands. The curve is the fit described in the text. 
The insert focuses on the low mass region. The dotted line in 
the insert indicates the D*j(2317) + mass. 
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FIG. 34: The invariant mass distributions of MC samples of 
(a) D sJ (2460) + and (b) £> s i(2536) + decay to D+tt+tt". The 
curves are the functions described in the text. 



• D sJ (2460)+ -*■ D+tt+tt- decay. 

• D sl (2536)+ -► .D+tt+tt- decay. 

• A third-order polynomial to describe the back- 
ground. 

Each signal decay mode is described by a PDF consisting 
of three Gaussians with a common mean, the parameters 
of which are determined using fits to MC samples (gener- 
ated with T = 0.1 McV). The fits for the D sJ (2460)+ and 
-D s i(2536)+ mesons are shown in Fig. 34 and correspond 
to mass resolutions of approximately 1.3 and 1.9 MeV/c 2 , 
respectively. The D s j(2460) + and Z? s i(2536)+ masses 
are allowed to vary in the fit and the D*j(2317)+ mass 
is fixed to the value determined in Section IX. 

The result of the likelihood fit is shown in Fig. 33. 
A L> sJ (2460)+ mass of (2459.7 ± 0.2) MeV/c 2 and 



TABLE VII: A summary of systematic uncertainties for the 
D s j(2460)+ and D s i(2536) + masses from the analysis of the 
DJ~Tr + Tr~ final state. 



Mass Uncertainty ( MeV/c 2 



Source 


D sJ (2460) + 


_D s i(2536) + 


Dj mass 


0.60 


0.60 


Signal line shape 


0.03 


0.01 


Background function 


0.01 


0.01 


Solenoid field 


0.05 


0.08 


Magnetization 


0.02 


< 0.01 


Energy-loss correction 


0.52 


0.30 


4> dependence 


0.01 


0.01 



Quadrature Sum 



0.80 



0.68 



D sl (2536)+ mass of (2534.3 ±0.3) MeV/c 2 are obtained 
from raw yields of 123 ± 15 and 193 ± 22, respectively 
(statistical errors only). A raw D*j(2317)+ yield of 6±3 
is also obtained, consistent with zero. 

The uncertainties in the D sJ (2460)+ and L> sl (2536)+ 
masses are summarized in Table VII. The largest uncer- 
tainty (0.6 MeV/c 2 ) is in the assumed D+ mass. In com- 
parison, uncertainties associated with the likelihood fit, 
such as the background shape or mass resolution, are rel- 
atively small. The remaining uncertainties are attributed 
to potential momentum biases in the tracking detectors. 

Uncertainties in tracking have several sources, dis- 
cussed in detail as part of a recent measurement of the A c 
mass from this collaboration [16]. A similar treatment of 
these uncertainties is reproduced for this analysis. The 
magnetic field of the BABAR detector is known to high 
precision, and variations of both the overall strength of 
the solcnoidal field and of magnetization effects of PEP-II 
beam elements produce no significant variation in mea- 
sured D+tt+tt - mass. Residual 0-dependent momentum 
tracking biases are also insignificant. The largest uncer- 
tainty arises from material inside the tracking volume 
and the effect of this material on energy loss corrections. 
Studies of large samples of A and Ks decays suggest that 
either the amount of material and/or atomic weight com- 
position of the SVT is underestimated by approximately 
20%. MC studies are used to estimate the bias introduced 
by this underestimation. A correction of 0.46 MeV/c 2 
and 0.29 MeV/c 2 for the L> sJ (2460)+ and L> s i(2536)+ 
masses is calculated based on these studies. The final 
results are: 



m(D s j(2460)+) 

to(L> s i(2536)+) 



(2460.2 ± 0.2 ± 0.8) MeV/c 2 
(2534.6 ±0.3 ±0.7) MeV/c 2 , 



where the first error is statistical and the second is sys- 
tematic. 

The D s j(2460)+ and L> s i(2536)+ signal PDF used in 
the likelihood fit of Fig. 19 was based on intrinsic widths 
of r = 0.1 MeV. After applying the same likelihood- 
integration technique described in Section IX for the 
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TABLE VIII: A summary of systematic uncertainties for the 
D s j(2460) + and -D s i(2536) + yields from the analysis of the 
D+ 



TABLE IX: A summary of the combined mass and width 
results. For the masses, the first quoted uncertainty is statis- 
tical and the second is systematic. The limits on the intrinsic 
width r are at 95% CL. 





Relative Yield (%) 












Particle 
Z?:, 7 (2317) + 
D sJ (2460)+ 
D s i(2536)+ 


Mass (MeV/c 2 ) 
2319.6 ± 0.2 ± 1.4 
2460.1 ± 0.2 ± 0.8 
2534.6 ± 0.3 ± 0.7 




Source 

Tracking Efficiency 
Signal line shape 
Background function 


LW(2460) + 
2.7 
2.7 
6.1 
3.9 
3.9 
9.1 


D s i(2536) + 
2.7 
0.3 
5.6 
3.9 
3.9 
8.3 


T ( MeV) 

< 3.8 

< 3.5 

< 2.5 


Df efficiency 
Df -> 07T+ yield 
Quadrature Sum 









D+7T final state, and after considering systematic uncer- 
tainties in the background shape and mass reconstruction 
resolution, the result is a 95% CL limit of V < 3.5 MeV 
and T < 2.5 MeV for the £> sJ (2460)+ and L> sl (2536)+ 
mesons. 

To estimate a total D*j(2317) + yield, a weight in- 
versely proportional to the p*-dependent D*j(2317) + se- 
lection efficiency is applied to each _D+7r + 7r~ combina- 
tion and the likelihood fit repeated. A similar process 
is used for the D sJ (2460) + and £> s i(2536)+ yields. The 
results are yields of 60 ±40, 760 ±110, and 1210±160for 
the D* sJ (2317)+, L> sJ (2460)+, and L> sl (2536)+ mesons, 
respectively, and for p* > 3.2 GeV/c (statistical errors 
only). 

The systematic uncertainties in the Z? s j(2460) + and 
D s i(2536) + yields are summarized in Table VIII. The 
cross section results are: 

a(D* sJ (2317)+ -► D+-K+-K-) = (0.3 ± 0.2 ± 0.0) fb 
a{D s .,{2A60) + -► D+tt+tt-) = (3.3 ± 0.5 ± 0.3) fb, 

S( J D s i(2536) + ^ J D+7r+7r-) = (5.2 ± 0.7 ± 0.4) fb, 

where the first error is statistical and the second is sys- 
tematic. 



XVI. COMBINED MASS AND WIDTH 
RESULTS 



TABLE X: A summary of yield results. All cross sections are 
calculated for p* > 3.2 GeV/c. The first quoted uncertainty 
for the central value is statistical and the second is systematic. 
The limits correspond to 95% CL. 



Decay Mode 



Central Value (fb) 



a(D* sJ {2317) + 

Dfn" 

Dfj 

D+n°ir° 

L>:(2112) + 

D+7T + 7T- 

cr(D sJ (2460) + 
Dfn" 
Dfy 

Dfn"-y [a] 
D*(2112) + 
D* sJ (2317) 
Dfn"n" 
D+77 
£>*(2112) + 

D+Tl + TT- 



)B(D* sJ (2317) 
115.8 ± 
-2.4 ± 
8.7 ± 

7.4 ± 
7 -0.5 ± 

0.3 ± 

)B(D sJ (2460) 

-1.0 ± 

14.4 ± 

42.7 ± 

7T° 41.6 ± 

1.1 ± 

5.5 ± 
3.5 ± 

-0.9 ± 
3.3 ± 



X)B(D. 



T~ 



Limit (fb) 



•) 



2.9 ± 
2.3 ± 

6.9 ± 
4.5 ± 
3.2 ± 
0.2 ± 



8.7 
8.9 
5.0 
2.2 
8.1 
0.0 



X)B{D+ -► 4>n+) 



7 



1.4 ± 

1.0 ± 

3.5 ± 

5.1 ± 
5.1 ± 

5.4 ± 
4.3 ± 

3.5 ± 
0.5 ± 



0.1 
1.4 

4.2 
5.0 
5.0 

2.4 
1.7 
4.1 
0.3 



< 15.7 

< 29.0 

< 20.6 

< 16.7 

< 0.6 

< 1.7 



< 15.2 

< 28.5 

< 13.2 
< 9.7 



a(Dl,(23l7)^)B(L 
Dfn+ 


>: J (2317) ++ - X)B{Dt - 


.07T+) 


< 1.5 


a(D:j(23l7Y')B(Dl 
D+ir- 


,(2317) u -» 


X)B(Dt -» 


dm 


+ ) 


<2.0 


a(D al (2536) + )B(D s 

Dfn+iv- 


l(2536)+ -» X)B(Df - 
5.2 ± 0.7 ± 0.4 


■^ + ) 


— 



[a] Includes both D*(2112) + 7r° and D* J (2317)+7 channels. 



The mass and width results for the D* 7 (2317) + , 
D s j(2460) + , and D s i(2536) + mesons are summarized in 
Table IX. The D*j(2317)+ and D s i(2536)+ mesons are 
observed in only one decay mode covered by this analysis; 
those portions of this table are copied unchanged from 
the respective sections of this paper. The D s j(2460) + 
mass is the average of that obtained from the Dfj, 
D+7r°7, and D+7r + 7r~ final states, although the latter 
measurement dominates in the average due to superior 
systematic uncertainties. This average is calculated us- 
ing a x 2 method that properly accounts for correlations 
among the systematic uncertainties in the three measure- 
ments. The limit on the intrinsic £> s j(2460) + width T is 
taken as the best limit obtained from these three decay 
modes. 



XVII. YIELDS AND BRANCHING RATIOS 

The eighteen decay yields a measured in this pa- 
per are collected in Table X. As described by Eq. 8, 
these numbers correspond to total yields for the decay of 
mesons having a center-of-mass momentum p* of at least 
3.2 GeV/c to a final state that includes a Df meson that 
decayed to 4>ir + . 

A 95% CL upper limit is calculated for those yields 
which are not statistically significant. These limits arc 
calculated using a frcqucntist approach by determining 
in each case the hypothetical value of a that is 1.96 stan- 
dard deviations above the measured values. In order to 
calculate the systematic uncertainty associated with a 
hypothetical value of a, those uncertainties that are pro- 
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TABLE XI: A summary of branching-ratio results. The first 
quoted uncertainty for the central value is statistical and the 
second is systematic. The limits correspond to 95% CL. 
For the hypothetical D* sJ {2317) ++ and D*j(2317)° mesons, 
an unknown additional factor from the ratio of production 
cross sections is involved. A lower limit is quoted for the 
D s ./(2460)+ -> D:(2112) + 7r results. 



Decay Mode 


Central Value 




Limit 


B(D2 j(2317) + -» 


X) /B(D* a j (2317)+ -D+7T U ) 
-0.02 ± 0.02 ± 0.08 


< 0.14 


D+7Y°7V 


0.08 ± 


0.06 ± 


0.04 


< 0.25 


D+rr 


0.06 ± 


0.04 ± 


0.02 


< 0.18 


L>*(2112) + 7 


0.00 ± 


0.03 ± 


0.07 


< 0.16 


D+TT + 1Y- 


0.0023 ± 


0.0013 ± 


0.0002 


< 0.0050 


B(D sJ (2460) + -> 
Dt*° 


X)/B(D a j (2460)+ — D+7T U 7 ) [a] 
-0.023 ± 0.032 ± 0.005 


< 0.042 


D+! 


0.337 ± 


0.036 ± 


0.038 


— 


D* s (2ll2) + TT a 


0.97 ± 


0.09 ± 


0.05 


> 0.75 


D: j (2317) + 7 


0.03 ± 


0.09 ± 


0.05 


< 0.25 


D + W°Tt° 


0.13 ± 


0.13 ± 


0.06 


< 0.68 


Df-fj 


0.08 ± 


0.10 ± 


0.04 


< 0.33 


D s *(2112)+ 7 


-0.02 ± 


0.08 ± 


0.10 


< 0.24 


DflV + TT- 


0.077 ± 


0.013 ± 


0.008 


— 


a{D* sJ {2317) ++ )/o{D* sJ {2317) 
B(D* sJ (23l7)++ -^ X)/B(D* sJ < 


+ )x 
;2317)+ -> 


DU°) 


< 0.017 


<j(D* sJ (2317) u )/<t 
B(D* S j(2317)° - 

D+7T- 


(D* sJ (2317) + )x 
X)/B(D* B j (2317)+ - D 


U°) 


< 0.013 



[a] Denominator includes both D*(2112)+7r° and D* J (2317)+7 
channels. 



portional to the signal yield (such as uncertainties related 
to selection efficiency) are scaled as appropriate. The re- 
sults are shown in Table X. 

The yields listed in Table X are used to calculate 
the branching ratios shown in Table XL For the 
Z?*j(2317) + meson, only one decay mode has been ob- 



served; this is used as the denominator when calculating 
the L>*j(2317)+ branching ratios. For the D sJ (2460) + 
meson, the Dfir°j decay mode (consisting of possible 



decay through either D*(2112)+7r° or D* sJ (2317) + j) is 
chosen for this role. For completeness, the yield from 
hypothetical D* sJ (2317) ++ -» D+tt+ and D* 7 (2317)° -> 
Dfir~ decays is compared to D* sJ {23l7) + — > Dfir° de- 
cay to produce a quantity that is proportional to both 
the respective branching ratios and production rates. 

In order to calculate the systematic uncertainty in the 
branching ratio results, systematic uncertainties common 
to the nominator and denominator are first discarded. 
Such common uncertainties include those associated with 
the Df — ► 4>ir + branching ratio correction and, in some 
cases, those associated with ir° and 7 reconstruction ef- 
ficiency. 

A 95% CL upper limit is calculated for those branching 
ratios associated with decay modes that lack statistical 
signficancc. The method used is similar to the frequcntist 
recipe used for the upper limits on the cross section. The 



results are included in Table XL 



XVIII. DISCUSSION 

The results in this paper confirm previous measure- 
ments, with generally higher precision. No new decay 
modes have been uncovered. Lacking any additional evi- 
dence to the contrary, the J p = + and J p = 1+ assign- 
ments for the D*j(2317) + and L> sJ (2460)+ mesons re- 
main a viable hypothesis. Except, perhaps, for the mass, 
there is currently no conflict with the interpretation of 
the L>*,(2317) + and L> s ./(2460) + mesons as the J p = 0+ 
and J = 1 + p-wave cs states, as shown in Fig. 1. 

The -D s j (2460) + meson mass measured in this analy- 
sis is an improvement over previous measurements [2-4] . 
High precision is obtained by using a decay mode that 
includes only charged particles coupled with a detailed 
understanding of the performance of the BABAR tracking 
detectors. The determination of the £>*j(2317) + mass 
remains limited by uncertainties in EMC energy scale. 

No intrinsic width of any statistical signficance is ob- 
served for cither the D* 7 (2317)+ or D sJ (2460)+ meson. 
There is sufficient detector resolution in the D + 7r + 7r~ 
decay mode that better limits for the D s j(2460) + meson 
should be attainable in the future with higher statistics. 
In contrast, limits on the D*j(2317) + width are currently 
limited by systematic uncertainties and will be difficult 
to improve. 

Both the mass and width results for the Z? s i(2536) + 
meson obtained from the decay to D + 7r + 7r~ are consis- 
tent with prior measurements [15] and theoretical expec- 
tations [5-8]. 

The analysis of the decay D s j (2460)+ -► £>*(2112)+7r° 
is complicated by the overlap in kinematics of the hy- 



pothetical decay D s j(2460) 



D; j (2317)+ 7 . The 



analysis presented here applies a likelihood technique 
to separate these two decay modes and provide a limit 
on the latter. A previous analysis from BELLE [3] 
of the decay £> sJ (2460) + -► D*(2112)+7r° did not ac- 
count for such possible contamination. As a consequence, 
it should be noted, their quoted yield for the decay 
£> sJ (2460) + -► D*(2112)+7r° should be compared to the 
D s j(2460) + — * Z? + 7r°7 result presented here in which 
both decay channels (L>* J (2317)+7 and D* s (2112)+tt q ) 
are included together. Specifically, the BELLE measure- 
ment of the branching fraction from e + e~ — > cc pro- 
duction is B(D sJ (2460)+ -► J D+7)/S( J D s . / (2460)+ -► 
L»*(2112)+7T ) = 0.55 ±0.13 (stat.) ± 0.08 (syst.) and 
can be compared directly to the result of this analysis 
(0.337 ± 0.036 (stat.) ± 0.038 (syst.)). The less precise 
BELLE result is somewhat larger, but still compatible 
within uncertainties. The result presented here is consis- 
tent with BELLE (0.38 ± 0.11 (stat.) ± 0.04 (syst.)) [12] 
and BABAR (0.275 ± 0.045 (stat.) ± 0.020 (syst.)) [13] 
measurements from D s ,/(2460) + production in B decay. 
The remaining branching ratio results shown in Ta- 
ble XI are generally competitive with the corresponding 
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BELLE limits [3] . An exception is the limit on the decay 
D* 7 (2317)+ -► L»+7. The BELLE publication, however, 
does not address the difficulties associated with modeling 
the D*j(2317)+ -> L>+tt reflection (see Fig. 14) which is 
the source of much of the systematic uncertainty quoted 
in this analysis. 

The searches for the D+7r°7r° and -D+77 decay modes 
of the D* 7 (2317)+ and L> sJ (2460)+ have been reported 
in this analysis for the first time. Both final states suf- 
fer from large backgrounds. More precise information on 
these modes will require more advanced analysis tech- 
niques, perhaps with the use of B meson decay or some 
other method of background suppression. 

The D* sJ (23l7) + and L> sJ (2460)+ masses are con- 
siderably lower than predictions from potential mod- 
els developed before their discovery. Since the 
-D* J (2317) + /D ;s j(2460) + mass splitting is so much larger 
than that for the -D s i(2536)+ and _D S 2(2573)+ mesons, 
it would appear that this conflict is intrinsic. Studies 
from Cahn and Jackson [17] suggest, however, that it 
is possible with some adjustment to fit the -D* 7 (2317)+ 
and Z? s j(2460)+ into a pcrturbative model. In any case, 
given the approximate nature of these models, there is a 
danger that even large discrepancies can be overstated. 

In contrast, according to Bardcen, Eichten, and 
Hill, models based on chiral symmetry naturally ac- 
commodate the large L>* 7 (2317)+/L> s j(2460)+ mass 
splitting [18] when the L>* 7 (2317)+ and D Si7 (2460)+ 
multiplct is treated as a chiral partner of the ground 
state D+ and D*(2112)+ mesons. Their model 
provides specific predictions for branching ratios, no- 
tably £(L> s j(2460) + -> Dfn + TT-)/B(D sJ (24:60) + -► 
D*(2112)+7r°) w 0.09 and £(D sJ (2460)+ -► 
D+j)/B(D s j (2460)+ -► D*(2112)+7r°) « 0.24, 
both in good agreement with our measurements of 
(0.077 ±0.015) and (0.337 ± 0.052), respectively. Other 
branching ratio predictions arc £>(.D*j(2317)+ — ► 
D*(2112)+j)/B(D* sJ (23l7)+ -► Dfn°) w 0.08 and 
BLD sJ (2460)+ - D s * i7 (2317)+ 7 )/B(D sJ (2460)+ -» 
D*(2112)+7r°) sw 0.13. Neither decay is observed by this 
analysis, but the corresponding upper limits still exceed 
the predicted branching ratios. This model also predicts 
a total intrinsic width of 23.2 keV and 38.2 keV for the 
.D*j(2317) + and D Si7 (2460)+ mesons, respectively, far 
below the limits measured here. 

Pcrturbative calculations from Godfrey [19] predict 
branching ratios of 0.19 and 0.55 for the radiative 
transition decay of the D* sJ (2317) + and D sJ (2460)+ 
to -D*(2112) + 7. These predictions are not consistent 
with the results reported here. This discrepancy could 
be considered evidence that both the D*j(2317)+ and 
D s j(2A60) + are not cs mesons. These theoretical pre- 
dictions, however, are plagued with large uncertainties in 
the partial widths of the isospin- violating _D* 7 (2317)+ — ► 
L>+7r° and L> s . 7 (2460)+ -» L>;(2112)+7r° decays. In ad- 
dition, leading-order corrections to the radiative tran- 
sition partial rates might explain current experimen- 
tal limits [20]. Nevertheless, the possibility that the 



-D*j(2317)+ andD s j(2460)+ are some type of new bound 
state, such as a four-quark DK molecule, should be taken 
seriously. 

The possibility of four-quark states has long been pro- 
posed [21, 22]. Perhaps the most unambiguous signature 
of the molecular interpretation of the -D* 7 (2317)+ meson 
would be the production of neutral and charged partners 
decaying to -D+7T . For an isospin 1 DK molecule, pro- 
duction of these isospin partners is expected at the same 
rate as the D* 7 (2317)+ [9]. This is clearly ruled out by 
our data. Other molecular interpretations, however, do 
not have this limitation [9, 10], and it may be challenging 
to rule them all out with the data from the i?-factories 
alone. 

Finally, now that more precise data on the -D* 7 (2317)+ 
and D s j(24:60) + is available, measurements of the other 
cs states have become more important. The measure- 
ment of the D s i(2536) + mass in this analysis will be 
helpful in further constraining models. With high statis- 
tics samples available from both B-factory experiments, 
improvements in the parameters of the Df, D*(2112)+, 
and D S 2(2573) + mesons should not be neglected and will 
hopefully follow soon. 



XIX. CONCLUSION 



An updated analysis of the D* sJ (2317)+ and 



cc 



LW(2460)+ mesons using 232 fb _1 of e+er 
data is performed. Established signals from the decay 
D*j(2317)+ -► D+tt° and L> sJ (2460)+ -► D+tt -/, L>+ 7 , 
and -D + 7r + 7r - are confirmed. A detailed analysis of in- 
variant mass distributions of these final states including 
consideration of the background introduced by reflections 
of other cs decays produces the following mass values: 

m(D*j(2317)+) = (2319.6 ± 0.2 ± 1.4) MeV/c 2 
m(D sJ (2460)+) = (2460.1 ±0.2 ±0.8) MeV/c 2 , 

where the first error is statistical and the second sys- 
tematic. Upper 95% CL limits of L < 3.8 MeV and 
L < 3.5 MeV are calculated for the intrinsic D* J (2317)+ 
and -D s j (2460) + widths. All results are consistent with 
previous measurements. 

The following final states are investigated: 
Dfir , D+7, i?: (2112)+ A J D; j (2317)+ 7 , D+7r°7r°, 
D*(2112)+7, D+77, D+7T+, and D+tt+tt". No sta- 
tistically significant evidence of new decay modes is 
observed. The following branching ratios arc measured: 



g(£> s , 7 (2460) + ->.D + 7) 

e(n 3j (2460) + -»_D + T°7) 

g(£> s j(2460) + ^£>+7r + 7r~) 

B(n 3 . 7 (2460) + -»n + ir 7) 



0.337 ±0.036 ±0.038 
0.077 ±0.013 ±0.008 



where the first error is statistical and the second sys- 
tematic. The data are consistent with the decay 
-D s j(2460)+ — ► D+7r°7 proceeding entirely through 
L»;(2112)+7r°. 



28 



Since the results presented here are consistent with 
J p = + and J p = 1 + spin-parity assignments for the 
D* sJ (2317) + and D sJ (2460) + mesons, these two states 
remain viable candidates for the lowest lying p-wave 
cs mesons. The lack of evidence for some radiative 
decays, in particular D* sJ (2317)+ -> D*(2112)+7 and 



D sJ (2460)- 



_D*(2112)+7, arc in contradiction with 



this hypothesis according to some calculations, but large 
theoretical uncertainties remain. No state near the 
D* sJ (231 7) + mass is observed decaying to Dfn ± . If 
charged or neutral partners to the D*j(2317) + exist (as 
would be expected if the D* sJ (2317) + is a four-quark 
state), some mechanism is required to suppress their pro- 
duction in e + e~ collisions. 
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